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ABSTRACT 

Many astrophysical binaries, from planets to black holes, exert strong torques on their 
circumbinary accretion disks, and are expected to significantly modify the disk struc- 
ture. Despite the several decade long history of the subject, the joint evolution of 
the binary + disk system has not been modeled with self-consistent assumptions for 
arbitrary mass ratios and accretion rates. Here we solve the coupled binary-disk evo- 
lution equations analytically in the strongly perturbed limit, treating the azimuthally- 
averaged angular momentum exchange between the disk and the binary and the mod- 
ifications to the density, scale-height, and viscosity self-consistently, including viscous 
and tidal heating, diffusion limited cooling, radiation pressure, and the orbital decay 
of the binary. We find a solution with a central cavity and a migration rate similar 
to those previously obtained for Type-II migration, applicable for large masses and 
binary separations, and near-equal mass ratios. However, we identify a distinct new 
regime, applicable at smaller separations and masses, and mass ratio in the range 
10 -3 < q < 0.1. For these systems, gas piles up outside the binary's orbit, but rather 
than creating a cavity, it continuously overflows as in a porous dam. The disk pro- 
file is intermediate between a weakly perturbed disk (producing Type-I migration) 
and a disk with a gap (with Type-II migration). However, the migration rate of the 
secondary is typically slower than both Type-I and Type-II rates. We term this new 
regime "Type-1.5" migration. 
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1 INTRODUCTION 

Understanding the co-evolution of binaries and 
cretion disks is fundamental in several fields of 
trophysics, including planet formatio n and migration 



(|Goldreich fc T rcmainc 1980; Ward 1993), patterns in plan- 
etary rings (|Goldreich fc Tremaind 119821 ). stellar binaries, 
compact object, and binaries involving supermassive black 
holes (SMBHs). 

Despite the long history of the subject, there are no 
self-consistent analytical models for the co-evolution of bi- 
naries and accretion disks, incorporating the fundamental 
physical effects over the long timescales on which the binary 
separation evolves. The standard q-mo del of radiatively ef- 
ficient turbulent thin accretion disks (jShakura fc Sunvaevl 
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Il973f ) relates the effective kinematic viscosity of the disk 
to the pressure v oc ap. The viscous evolution of the disk, 
however, is often modeled wi thout considering the pressur e 
dependence of the viscosity (|Lvnden-Bell fc Pringlel[l974l ). 
Similarly, models of the gravitational interaction between 
the disk, which describe the launching of spiral density waves 
in the disk that remove angular momentum from the binary, 
also do not account for the tidal heating of the disk and the 
corresponding feedback on t he torque cutoff phenomenon 
(|Goldreich fc Tremainelll980r ). 

The evolution of the circumbinary disk is sensitive to 
the above mentioned assumptions, especially when the mass 
of the secondary is large, and can strongly perturb the disk. 
For a massive secondary, the tidal torque clears a gap in 
the disk, and the viscous radial inflow of the gas pushes 
the object inward on the viscous timescale (Type-II migra- 
tion). If the secondary mass, m B is larger than the local disk 
mass, rrid = 47rr- 2 E, where E is the surface density, then the 
migration slows down, as the spiral density waves cannot 
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remove angular momentum away from the binary at a rate 
on which the gas flows in. This leads to the pile-up of gas 
outside the secondary's orbit, in which the gas density in- 
creases by up to a factor £>~ 3//8 , where B — mdo/m B < 1 
and r wdo is the unperturbed local disk mass (|Sver fc Clarke! 
Il995h . Once this steady-state level is reached, the viscous gas 
inflow velocity matches the inward migration of the object 
(secondary dominated Type-II migration). 

In this paper, we foc us on such systems, w ith m s > m^o, 
and point out that the ISver fc Clarke! l| 19951 ) steady-state 
level of gas pile-up cannot be reached for sufficiently large 
secondary masses ?n s 3> m d o, where B < 1. The en- 



-5/8-, 



can increase 



hanced viscous dissipation rate (D„ oc B 
the disk temperature such that it becomes radiation pres- 
sure dominated. The enhanced pressure makes the disk puff 
up (H oc B~ 5 ^ s ), and reduces the relative gap size. Once 
the gas approaches within a distance less than the scale- 
height from the secondary, t he torque that the disk ex- 
erts on the binary h a s a cutoff (|Goldreich fc Tremaine|[l980l ; 
lArtvmowic3ll993bl lal; lGoodman fc Rafikovll200ll ) which lim- 



its the migration rate of the secondary. Once the gas en- 
ters the Hill radius, it can furthermore flow across the sec- 
ondary's orbit along horse shoe orbits or accrete onto the 
secondary. We derive an analytical quasi-steady-state model 
for the co-evolution of the disk and the orbital migration of 
the se condary, in which we combine a IShakura fc Sunvaevl 
(| 19731) disk with the theory of the binary-disk interaction 
bv lGoldreich fc Tremaine] (|l98fj ) self-consistently. In partic- 
ular, we adopt the viscosity prescription of standard thin ac- 
cretion disks proportional to pressure^] calculate the sound 
speed and vertical balance including both gas and radi- 
ation pressure (p gas and Prad), adopt the simple analyti- 
cal approximation to the angular m omentum exchange be- 
tween the binary and the disk of lArmitage fc NataraianI 
|2002l ). c onsider the standar d viscous and tidal heating of 
the disk (|Lodato et al.ll2009l ). and self-consistently account 
for the feedback on the pressure, viscosity, scale-height, 
and the torque cutoff near the s econdary's orbit. We gen - 
eralize the steady-state m odel of iHourigan fc Ward! l| 19841 ). 
IWard fc Houriganl j 19891 ). and lLiu fc Shapirol (|2010l ) by self- 
consistently including variations in the viscosity and pres- 
sure caused by the pile-up. We derive azimuthally aver- 
age d steady-state a n alytic al disk models which recover 
the iGoodman fc Tar] (|2004l ) solution for arbitrary /? = 
Pgas/(p g as +Prad) in the limit that the secondary mass m B 
approaches zero, but the disk structure is significantly mod- 
ified by the secondary over multiple accretion timescales for 
larger m s . 

The disk structure in this overflowing state with a 
pile-up is intermediate between the weakly perturbed case 
without a secondary and the case with a gap. Not sur- 
prisingly, the migration rate in such an intermediate state, 
which we label Type-1.5, is significantly different from the 
corresponding limiting cases of Type-I migration and the 
secondary-dominated Type-II migration. The transition be- 
tween Type-I and Type-II migration as a function of the 



1 Here v oc p ga s and v <x (p ga s + Prad) are respectively known 
as a and /3-models. We formulate our problem for a general a 
or /3-disk in Sec. 2, but then derive the analytical results for the 
special case of a /3 disk viscosity. 




Figure 1. Gas pile- up and overflow in a circumbinary accretion 
disk with component masses M. and m B , binary separation r B , 
and accretion rate M. We distinguish five distinct radial zones: 
an inner and an outer far zone where the effects of the secondary 
are negligible, an interior and an exterior near zone where the 
tidal effects are significant, and an extended middle zone with a 
significant gas pile- up (see §H). 



secondary mass was previously typically investigated by 
considering only the change in the surface density due to 
gap formation, but without investigating the feedback from 
the changes in viscosity and pressure (|Hourigan fc Wan 



1984; 



Ward fc Houriganl Il98 
Ward! 1 19971 : iBate et ail 



ixd 
ou 



Korvcanskv fc Papaloizo 
20031 : ICrida fc Morbidellil 
simulations show that migration is sen- 
sitive to temperature variations and radiation pressure 
(iD'Angelo et al.ll2003l : iPaardekooper fc Mellemall2006l . 120081 : 



IKlev fc Cridall2008l ). We derive the Type-1.5 migration rate 
for the self-consistent radial profile including these effects 
when the pile-up is significant in an overflowing steady-state 
disk. As the secondary migrates inwards across the increas- 
ingly hotter inner regions of the disk, the gap opening con- 
ditions and the migration rate change even if one neglects 
the feedback on viscosity and temperatu re due to gas pile-up 
llHaiman et al ] |2009l :lK ocsis et aDl201ll ). but here we show 
that the changes are significant over a much wider range of 
masses and radii in the self-consistent model. We discuss mi- 
gr ation and gap open ing for SMBH binaries in more detail 
m iKocsis et aTJ (|2012l ), hereafter Paper II. 

The remainder of the paper is organized as follows. In 
§ [2] we lay out the basic equations governing the hydrody- 
namical and thermal evolution of the disk, as well as the 
migration of the secondary. We solve the equations numer- 
ically in § [3] We then derive an analytical solution in § [4] 
We summarize the results, and discuss how they depend on 
the most important physical parameters, in §[5] We offer our 
conclusions in §|6] A more detailed discussion and the impli- 
cations for SMBH binary systems is presented in Paper II. 

We use geometrical units G = c = 1, and suppress fac- 
tors of G/c 2 and G/c 3 to convert between mass, length, and 
time units. Our basic notation for the disk and secondary 
parameters are depicted in Figure Q] 
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2 THERMO-HYDRODYNAMICAL 

INTERACTION BETWEEN A DISK AND A 
SECONDARY 

We examine the evolution of the secondary and an az- 
imuthally and vertically averaged Shakura-Sunyaev disk (i.e. 
axisymmetric one-zone disk) in local thermal equilibrium. 
Here we review the basic equations. First, we write down the 
continuity and angular momentum transport equations in- 
cluding the viscous torque and the gravitational tidal torque 
of the secondary. The back-reaction of the tidal torque 
changes the angular momentum of the secondary. The vis- 
cous and tidal torques depend on the disk surface density, 
viscosity, and pressure gradient (or scale- height). We derive 
the vertically averaged disk structure assuming that: (i) the 
local viscous plus tidal heating equals the radiative cool- 
ing with photon diffusion limited vertical radiative flux (i.e. 
negligible radial heat transport); (ii) the viscosity is propor- 
tional to either gas+radiation pressure (a-disk) or just the 
gas pressure (/3-disk), and (iii) gas plus radiation pressure 
supports the disk against the vertical gravity. This yields a 
closed set of nonlinear partial differential equations for the 
disk and the location of the secondary in 1+1 dimensions. 
We present solutions in subsequent sections below. 



2.1 Angular momentum transport 

We denote the masses of the primary and the secondary 
objects by M, and m s , the surface density of the disk by E 
(assuming axisymmetry) , and the radial bulk velocity of the 
disk by v r , which is negative if gas accretes toward r = 0. 
The continuity and angular momentum equations for the 
disk arfl 

= 2nr <9 t E + d r (2nrEv r ) , (1) 
d r T = 2-Kr 9i(Er 2 f!) + d r (2Trrv r 'Er 2 Q, ) , (2) 

where the total torque T — —T v + Td is due to viscosity and 
the gravity of the secondary, given by 

T v = -27rr 3 (<9 r fi)^E ~ 3irr 2 flvZ, (3) 



d r T d = 27rrAE . 



(4) 



Here A is the torque per unit mass in the disk, approximately 
given by 



A 



where 



-i/gVfW 4 /A 4 if r < r s 
+ i/g 2 r 2 tt 2 r 4 /A 4 if r>r s 



A = max(]r — r s \, H) 



(5) 



(6) 



q = m s /M m , H <C r is the scale-height of the disk, and 
/ is a constant calibrated with sim ulations. This approxi- 
mate formula for A, introduced by lArmitage fc NataraiarJ 
(|2002l ). accounts for the net contribution of all Lind- 
blad resonances as well as the torque cutoff within r s ± 



2 In our notation, S r S = 8T,/dr and dtT, = 8T,/dt = E. T 
refers to torque, T v and T d are viscous and tidal torques as in 



IChang et ail j201dh : IlIu fc Shapirol d201Ch . The angular momen- 
tum flux is Fj = T. Central and surface temperatures are labelled 
with T c and T s . 



H (|Goldreich fc Tremaind Il980l ; IWardl 1 19971 ). and guar- 
antees that the torque vanishes at r S> r s . Here / = 
(32/8l7r )[2jf (2/3) + K 1 (2/3)] 2 = 80 outside the torque 
cutoff in Goldreich fc Tremaind (|l980T ). / = 0.23 x (3/2ir) = 
0.11 in iLin fc Papaloizoul(|l986E . an d / = 10 2 calibrated to 
match the gap opening conditions in lArmitage fc Nataraianl 



(200 



We adopt a conservative value /_ 



f/w- 



in our numerical calculations, but keep the /_2 terms gen- 
eral in all of our analytical formulas. Note that practi- 
cally Eq. (fSj assumes that the tidal torque density "sat- 
urates" instead of having a true cutoff near the sec- 
ondary as long as t he gas density is non-vanishing there 
|Artvmowiczl Il993bh . which accounts for the effects of 
shocks nea r the secondary ((Go odman & Rafikov 1200 ll ; 
iDong et all l201ll ; iDuffell fc MacFadvenl |2012D I 4 I l 5 l How- 
ever, this prescription might be inaccurate for a high- 
mass secondary forming a gap in the disk, where the tidal 
torques are due to spiral str eams passing near the sec 



onda r y on horse-shoe orbits dMacFadven fc Milosavlievic 
20081; IShi et al.l [20121; iRoediget al.l |2012l: iBaruteau et al 



20121 ; iPetrovich fcR afikov 20l2). We do not consider the 



torques inside the Hill radius, \r — r&\ < j"h = (<z/3) ' r s , as- 
suming that gas reaching this region flows in across the sec- 
ondary's orbit. Outside this region, we use Eq. Q, assume 
that gravity is dominated by M. , and the orbital velocity is 
nearly Keplerian, Q ~ M.f~ 3 ^ 2 , where f = r/M,. 

After some algebra (|Frank et al. I l2002tl . Eqs. (HHSJ) sim- 
plify to 



27rr 



Or 



d r T 



d r (r 2 D.) 



d r T 



2nrY,d r (r 2 n) 



(7) 



The total mass flux across a ring of radius r is defined as 



M(r, t) 



-2TTr'Sv r 



d r T 



(8) 



d r (r 2 D.) ■ 

Eq. along with the definition of the total torque T in 
Eqs. ((3}{5j) describes the evolution of the axisymmetric disk 
surface density and radial velocity as a function of radius 
and time. 

The evolution of the secondary's orbital radius, r s , is 
driven by the tidal torques of the gas and gravitational wave 
(GW) losses. The angular momentum of the secondary is 
L s — m B r 2 fl s so that 



d r T d dr - T G 



(0) 



3 iLiu fc Shapirol ( |2010h used Eq. J5} with / = lQ- 2 . IChane et al.l 
ll2010h adopted a to r que m odel, extrapolating Eq. (18) of 
iGoldreich fc Tremaine! l ll980l ) (with a modified constant prefac- 
tor of / = 0.1 X (4/9-7T) ~ 10 -2 ), such that their torque den- 
sity approaches a constant at f > r s Th e linear perturba- 
tive analysis of IGoldreich fc Tr emaind dl98Cf) is not applicable 
if q > a 2 or q > 



ll997f[Korvcanskv fc Papaloizoulll99e ) 



Idrcicn & Ircniamc ( 198U ) is not applicable 
H/r) 3 llMever-Vernet fc Sicardvlll987l; IWardl 



Recent simulations llDong et all l201ll; iRafikov fc Petrovichl 
2011 IDuffell fc MacFadvenl l2012F have shown that the original 



Goldreich fc Tremaind lll980l ~ torque density is correct close to 



the secondary, but the actual torque decreases in amplitude and 
changes sign outside of r B + 3H. However the relative contribution 
of these outer regions to the total torque is negligible. 
5 We do not account for relativistic corrections to the tidal torque 
which are expected t o be small at t he separations in the gas-driven 
regime r s > 100M. l lHiratall2011al lbl). 



© 0000 RAS, MNRAS 000, 000-000 



4 Kocsis, Haiman, & Loeb 



where — Td is the recoil due to the torque exerted on the 
disk, Eq. @, and with f B = r s /M m , the torque from the 
GWs is given by 



2g\ 



32 mi . 



-7/2 



(10) 



Given v(r,t) and H(r,t), Eqs. {7} and ([9]) provide three 
equations for the three unknowns: £(r, t), v r (r, t), and v S i(t). 

We examine steady-state solutions to these equations 
where E = and dM/dr = so that M(r,t) = M is 
a constant[f] Note that in general the disk need not be in 
steady-state. However, in many cases the inflow rate of gas 
may be much faster than the radial migration speed of the 
secondary \v sr \ <C \v r \. If this is satisfied in a wide range of 
radii up to the outer edge of the disk, then, the secondary 
is effectively stationary in the azimuthally averaged picture, 
and the radial profile of the disk might be expected to re- 
lax to a steady-state, independent of the initial condition 
of the disk. We propose that the secondary then migrates 
slowly through a sequence of quasi-steady-state configura- 
tions of the disk with a fixed M(r, i) — const. Then, Eq. ([8]) 
becomes 



d r T v - d r T d = Md r (r 2 Q) . 



(11) 



This is a first-order ordinary differential equation for T v (r), 
once d r T d (r) is specified for a specific disk model. 

2.2 Boundary conditions 

We distinguish two types of inner boundary conditions cor- 
responding to whether or not the perturbation is strong 
enough to lead to a truncated disk with a wide hollow cir- 
cular cavity. Here, "wide" means wider than the Hill radius 
(see below). 

I. If E(r) all the way to the innermost stable circular 
orbit nsco of M. (i.e. the disk does not h ave a cavity) , we re- 
quire a zero-torque boundary condition (INovikov fc Thornej 
1 19731 ; IPenna et al]|2010l ; iTanakallioTH ; IZhu et al.ll2012i ). 



r^(nsco) = . 



(12) 



Starting with this boundary condition, we obtain, among 
other properties of the steady-state disk, the gas velocity 
profile v r (r). As stated above, if this inflow velocity is much 
faster than the migration velocity of the secondary over a 
large range of radii, then one might expect that the disk ap- 
proaches this steady-state configuration, independent of the 
initial condition. In the opposite case, the steady-state as- 
sumption may be violated by the time-dependent migration 
of the secondary. As we will show, the steady-state solution 
with a fixed M requires a large build-up of gas outside the 
secondary for the viscosity to overcome the tidal barrier of 
the secondary. We refer to these solutions, in which the disk 
is not truncated outside the secondary as "overflowing" . 

II. If the tidal torques dominate over the viscous torques 
near the secondary, gas is expelled from the region near the 
secondary and a wide gap forms. Assuming that the char- 
acteristic radius r g , where the tidal torque is exerted on the 
disk near the edge of the gap, tracks the inward migration 



6 As stated above, we negl ect the accretion o nto the secondary 
for simplicity (however, see lLubow et al1ll999lV 



of the secondary with r g = Ar s where 1 < A < 3 is a con- 
stant, we require that the gas velocity a t that radius satisfy 
(|Sver fc Clarkdll995l ; llvanov et al]|l999T ) 



v r (\r s ) = A?; s 



(13) 



Note that A is not specified by hand ab-initio; it is found by 
assuming steady-state in our solutions below. This condition 
can be understood intuitively, since the secondary cannot 
"run away" and leave the outer disk behind (if it did, it 
would cease to be able to torque the disk and would have to 
slow down). Likewise, the gap edge cannot get closer to the 
secondary (if it did, gas would pile-up and the gap would 
eventually close). Although the disk is not in steady-state 
near its boundary, we assume M(r) ~ M at r > r g (see 
discussion in § l4,1.4l below'). 

Based on Eq. (J5J) and ©, Eq. JT3J) is equivalent to 



d T d dr = m * r ° Q sM 
47rr2E(r g ) 



(14) 



Note that here and throughout the paper by "gap" 
we refer to situations where the gas density becomes ef- 
fectively zero outside the secondary, such that the inflow of 
gas from the outside pushes the secondary inward accord- 
ing to In these cases, we assume that inflow across the 
orbit is insignificant, and in particular, we neglect torques 
from the gas interior to the orbit. In our calculations, a gap 
is effectively a hollow circular cavity in the disk, which is 
supported by the tidal torques of the secondary. However 
we emphasize that we do not rule out the presence of a lo- 
cal density decrement, resembling an annular gap, with a 
significant mass flux across the gap. 

In practice, we attempt to find a solution with either of 
the above two boundary conditions, and then check whether 
the solution is self-consistent. By construction, only one of 
the two boundary conditions will lead to a self-consistent 
solution as confirmed below. 



2.3 Physical conditions in the disk 

Next we derive H(r) and u(r) which appear in the tidal and 
viscous torques in Eqs. ([3][5]). 



2.3.1 Vertical balance 

Let us first derive the scale-height, H. If the vertical gravity 
is dominated by M., (i.e. \r — r s \ > t~h), then in vertical hy- 
drostatic equilibrium H — c s /fi where c s — \fp[p is the lo- 
cal midplane sound speed andp = p gaa + p ra d is the pressure 
due to the gas and radiatiorjj, p gas = pkT c /(fj,m p ), p ra d = 
|aT 4 , where T c is the central temperature, a — 4a /c is the 
radiation constant, a is the Stefan-Boltzmann constant, m p 
is the proton mass, and /j, — 0.615 is the mean particle mass 
in units of m p . Since p = E/(2H), c\ = 2#p/£ = 2c s p/(En) 
so that c s = 2p/(EO). The pressure can be expressed as 



7 Note that the gas is not degenerate and is not isentropic, 
therefore the assumption of p oc p 5 / 3 or p 4 / 3 made in most nu- 
merical simulations of accretion disks is inappropriate. In fact, 
p oc T 4 oc M 4 p e for a radiation-pressure dominated standard 
Shakura-Sunyaev disk with no secondary. 



© 0000 RAS, MNRAS 000, 000-000 



Gas pile-up, gap overflow, and Type 1.5 migration: general theory 5 



P — Prad/(1 — P), where ft = p sas /p. If photons are trans- 
ported to the surface by diffusion then the mean radiation 
flux is 



3 



Here T 3 is the surface temperature, r = kE/2 is the 
optical depth from the midplane to the surface, where 
k — 0.35 cm 2 /g is the opacity assumed to be domi- 
nated by electron-scattering. We do not investigate changes 
caused by free-free opacity at large radii for simplicity. 
and neglect deviations from blackbody radiation (see e.g. 
iTanaka fc Menoul |20ld for more a detailed model). Thus, 
Prad = |kS_F/c, so that c s = nc^ 1 F Q^ 1 (1 — /3) _1 , and we 
have 



H= Sl = _H L 

n en 2 1 - 



(16) 



Note that Eq. (|16|l is valid in general for radiation flux lim- 
ited, geometrically thin disks, independent of the source of 
dissipation and viscosity|f| 



2.3.2 Viscosity 

In the standard Shakura-Sunyaev a and /3-disk models, the 
viscous stress tensor, ty = pvViVj satisfies t r<t> = — |a/3 p, 
whe re b — or 1, respectively, and a is a constant parame - 
ter (|Shakura fc Sunvaevlli973l ; ISakimoto fc Coronitilll98ll) . 
implying that 



v = ac s Hf3 = a- 



/3 b F 2 



(17) 



c 2 D. 3 (1 - /3) 2 ' 
In the second equality, we have substituted Eq. (|16l) . 

2.3.3 Local thermal equilibrium 

We assume steady-state thermal equilibrium in which heat 
generated by viscosity and the dissipation of the spiral den- 
sity wave escapes the optically thick disk in the vertical di- 
rection by photon diffusion. The vertical radiation flux is 
F = D v + Dd. The viscous dissipation rate per disk face 
element is 



A, = ^^ = V,E. 
47rr 8 



(18) 



We assume that the density waves generated by the tidal 
torque are dissipated locally in the disk and turned into 
heat, yielding the rate Dd- This is expected to be an 
adequate approximation ba sed on analytical arguments 
l|Goldreich fc Tremaind Il980l. Eq. 97 therein) and numer- 
ical studies jPong et all |201 ll ; iRafikov fc Petrovichl |2012| ; 
iDuffell fc MacFadvenlbOia T especially in the regime where 
the disk is strongly perturbed. 



8 One possible source of inconsistency is that convective verti- 
cal heat transport is conventionally neglected here. This may 
be significant for optically very thick, radiation pressure dom- 
inated disks (especiall y so-called /3 disk s) with a large vertical 
temperature gradient llBlaes et alj|201ll) . The heat transport in 
this regime may be analogous to the convection zones of stars. 



Following Goodman fc Rafikovl (|200ll ) and 
lLodato et alT(|2009l 1 9 l. 

Cd= (Q-Q)»r d = i 

47rr 2 
The total vertical flux or total dissipation rate is 

F = D v + D d = -Q, 2 vY> + i(n B - fi)AE. (20) 
8 2 

Using the above equations we derive E and T c for a given 
D v and F at each radius (see Appendix |A]) . 



2.3.4 Summary 

Combining the previous expressions, we obtain 
v _ 8 K/ 'ft) 4 ' 5 a 1 ' 5 p 1 -"^ Dt /5 



39/5 a 4/5 K l/5 



^4/5 pl/S 



(21) 



(Hmv/k) 1 ' 5 ^ 5 P (1 -" )/5 F x/s D i/s ^ (22) 



(23) 



31/5 Q l/5 a l/5 Ql/5 

where 

/3 (1/2)+(b - 1)/lu = c[fc/( At m p )] 2 / 5 9/10 Dl /W 

1-/3 (3ao-)Vio K 9/io 779/10 

All other disk parameters can be derived from these rela- 
tions. For example, the scale-height H and the quantity i/E 
that determine the torque (Eq. Hip are given by Eqs. (|16|l 
and (|18|l . In particular, the limiting cases for H are 



KC _1 n- 2 F 
y/k/(jimp)nT c 



1/2 



If P < 1 ; 

if $ ~ 1 . 



(24) 



In the limit that the only source of heat is viscosity in a Ke- 
pl erian disk, J 1 = F) v = (3/8ir)A'IQ , we recover the solution 
of iGoodmanl (|2003l ) up to a constant of order unitvF"! 

More generally, Eqs. (fl6]i. (pT|) . and J23J, along with 
the definition of D v and F in Eqs. (jTHJ and and the 

angular momentum flow equation (|11|) provide a closed set 
of equations for the stationary disk, valid throughout the gas 
and radiation-pressure dominated regions for a and /3 disks. 
The solution is self-consistent if for all r, the disk is thin 
(H < r), the radiation flux is sub-Eddington (L ~ 2-rrr 2 F < 
A-rcGM, / k), the radial accretion velocity is subsonic (v r — 
M/2-7rrE < c s = HQ), radial heat transport is negligible, the 
self-gravity of the disk is negligible and the disk is stable 
against fragmentation (Q — c 3 Q/ (nCE) 1), the disk is 
optically thick (r = hE/2 ^1), and the boundary conditions 
are satisfied (implying in particular that v al <C v r across a 
wide range of radii for overflowing solutions, see 



9 We add a factor of 2 that appears to be missing in lLodato et al.l 

this enters because of the two disk faces. 

10 We find a small difference in the dens i ty and temperature nor- 
malization constants, due to IGoodmanl l l2003f ) neglecting a 4/3 
prefactor in the vertical diffusion equation F = ^aT^/r. 

11 The model is furthermore self-consistent only outside the 
secondary's Hill sphere since the gravity of the secondary is ac- 
counted for as a perturbation to the primary's gravitational field, 
and the equations are linearized in the derivation of the torque 
formula. The tidal torque model is n evertheless often interpolated 
to within this region, as well (e.g. iGoldreich fc Tremainelll980l: 
lArmitage fc Nataraia r] |2002h . Here we avoid this extrapolation 
by excising the region within the Hill radius from our domain, 
assuming that gas entering this region flows across the secondary 
orbit. 
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We verify that these conditions are indeed satisfied for the 
overflowing solutions below. 



to the Hill radius as done p revio usly in Lin fc Papaloizoul 
<| 1986ft . ISver fc Clarke! (| 1995ft . and lLodato et all (|2009ft .~ 



DISK STRUCTURE 
SOLUTIONS 



NUMERICAL 



First we generate numerical steady-state solutions for tidally 
and viscously heated disks assuming that the migration rate 
is much smaller than the radial accretion velocity in the disk. 
These numerical solutions are useful to verify the detailed 
analytical estimates presented in the following section. 
We proceed along the following steps: 

(i) Obtain the ratio of gas to total pressure, /? = 
/3(r, D v , F), by inverting Eq. (|23p . A unique solution is guar- 
anteed by the intermediate value theorem, since the left 
hand side is a monotonic function of /3, mapping < /3 < 1 
to all positive real numbers, while the right hand side is 
positive and independent of ft. 

(ii) Substitute the solution for /3 in Eq. (fl6|) and Eq. (pTj) 
to obtain H(r, D V ,F) and E(r, D„,F). 

(iii) Substitute /3, H, and E in the definition of F, Eq. PCT)) 
to get an equation between F and D v for fixed r and r s . 
Invert this relation to find F(r,r B , D v ). Similar to step (i) 
one can show that the solution exists and is unique. 

(iv) Using Eqs. fTHHISJ, obtain the function d r Td = 
9d(r, r B ,T u ). 

(v) Substitute into Eq. 1)111) . to obtain an expression 
d r T v = g v (r,r B ,T v ) for a fixed M. Solve this differential 
equation for T„(r, r s ). 

(vi) Substituting back into D„ and F, Eqs. (fig)) and (f2U|) 
and the formulas of step (ii) to get E(r, r s ), T c (r,r s ), and 
H(r,r s ). 



(iii; 



m steps 



(i) 



The complexity is related to the nonlinearities in steps (i 
and|(v)| Nevertheless, the solution exists and is unique 
and 



(hi; 



(v) is a boundary 



However, step 

value problem of a nonlinear first-order differential equa- 
tion, which can have many solutions. We solve the differ- 
ential equation numerically upstream from an initial value 
T v (risco) = 0. Without the secondary, the solution is simply 
~ r ?sco^isco), which leads to the Shakura- 
Sunyaev disk. If q <C 1, then the secondary creates a small 
dip in T v (r) in its neighborhood, where the depth of the 
minimum increases with q. For larger q, T v (r) becomes very 
small positive approaching the secondary from downstream, 
and the surface density approaches zero. In this regime, tidal 
heating dominates over viscous heating, and H > \r — r s \, 
implying that the pressure gradient shifts the torques out of 
resonance, and the torque is suppressed according to Eq. ©. 
Since the adopted torque model is valid only outside the sec- 
ondary's Hill radius, we stop the calculation at r s — rn, and 
restart it at n = r s + ru assuming thaQ 

T v (r s - rn) « T„(n) . (25) 

This has a similar effect to smoothing the torque interior 

12 The tidal torque is monotonically increasing and decreasing, 
interior and exterior the secondary, respectively. The solution is 
uniquely determined by the initial value T^(risco) and T v (r r mi) 
in the two domains. However, r; can be arbitrary as long as rg > 



The solution is approximately self-consistent if the mi- 
gration rate is slower than the radial gas velocity outside 
the secondary. However, if this is not satisfied, a cavity 
opens and the disk becomes truncated. In this case, we seek 
a different solution in step |(v)| which satisfies the boundary 
condition in Eq. (|13l) . This is possible by increasing n in 
Eq. (gSJ where T v (n) ss 0, until Eq. J37|) is satisfied. Here 
ri can be identified as the truncation radius at the inner 
edge of the disk. We distinguish the characteristic trunca- 
tion or gap radius to reside at r s where the tidal effect is 
exerted on the disk, more specifically the boundary where 
the tidal torque density becomes subdominant and use r g 
in the boundary condition, Eq. (|37|l The surface density 
increases rapidly within n < r < r g has a maximum and 
decreases thereafter. We assume that the disk is truncated 
interior to n if a gap forms with n > r s + ru . 



DISK STRUCTURE 
SOLUTIONS 



ANALYTICAL 



Here we derive an analytical solution to the nonlinear equa- 
tions in § [2] Such solutions can be derived asymptotically 
far from the secondary or near the secondary, where either 
the tidal torque or the viscous torque dominates, or where 
the angular momentum flux is negligible. We therefore dis- 
tinguish the corresponding far, middle, and near zones (see 
Figure [TJ . The far zones are well inside and well outside the 
secondary, where the effects of the secondary are negligible. 
The middle zone is the region outside the secondary where 
the tidal effects (i.e. torque and heating) are locally neg- 
ligible compared to the viscous effects, but where the gas 
pile-up is significant and the disk profile is modified. The 
near zones are just inside and just outside the secondary's 
orbit, where the tidal effects of the secondary dominate over 
the viscous effects. We restrict the near zone to outside the 
Hill radius, where the adopted tidal torque formula is valid. 
In addition to providing a basic understanding of the disk 
structure, the approximate analytical solutions allow us to 
infer the migration rate of the secondary. 

To keep track of the approximations and notations in- 
troduced for the various zones below, we provide Table [T] 
for convenience. Note that the far/middle/near zones divide 
the disk into five radial slices, and the asymptotic behav- 
ior further depends on whether the disk becomes truncated 
forming a wide gap (in the middle zone) and whether the 
torque is saturated by the condition on the radial distance 
from the secondary is 5r = r — r s < H (in the outer near- 
zone). Each row in the Table corresponds to one of these 
disk regimes, discussed in a corresponding subsection below, 



13 In practice, we generate solutions for many different r\. We 
seek the radius r g at which the tidal torque cuts off in the numer- 
ical solution: T d (r g ) = 0.1 T d (r poak ) where r g > r pcak and r pcak 
is where T d (r) attains its maximum. We use this value as the 
gas velocity iy(r g ) in Eq. 137> . We find that the gas velocity is 
nearly constant in the neighborhood of r g and the surface density 
is near its peak, so the solution is insensitive to the details of this 
convention. 
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Md r (r 2 n) 


T„ 


d r T d 


Far zone 




/ 


<p(r, risco) 





Mid. with gap 


|4.1.3| 





T^(r B ,r e ) 





Mid. overflow 


|4.1.2| 





T b m c > s ,n) 





Near ext. uns. 


|4.2.2| 





CO, n,n) 




Near ext. sat. 


|4.2.3| 





ip{r, r B , n) 


rf/tf 4 


Near interior 


14.2.11 


/ 





-r 4 /|<5r| 4 



Table 1. Approximations and notations for the various radial 
radial zones in the disk, used in § [4] 



and shows which terms are relevant in Eq. The sub- 

dominant terms are marked with a "0" . The column with T v 
shows functions we introduced related to the viscous torque, 
and d r Td shows the scaling of the specific tidal torque in 
Eq. ©. 

In the following we mostly focus on /3-disks (i.e. 6=1) 
and examine both radiation and gas pressure dominated 
disks, but it is straightforward to derive analogous formu- 
las for a-disks in the same way. We also note that in the 
radiation-pressure dominated regime, the viscosity of a- 
disks is larger by a factor of p gas /(pgas + Prad) = P~~ ■ This 
would generally lead to stronger overflows for a smaller gas 
pile-up, and the cavity would close for a wider range of 
parameters than we find below for /3 disks. 



4.1 Far and middle zones 

First we examine the region sufficiently far from the sec- 
ondary, either inside or outside of its orbit, where 



\d r T d \ < <9 r T„ « Md r (r 2 n) . 



(26) 



In this region, Eq. (|lip can be integrated and substituted in 

mm 



T v = Mr Q + T bc 
3 fi 



F&D V 



8tt r 2 ' Tv 8tt 



MQ 2 + T bc 



(27) 
(28) 



where T^ c is an integration constant determined by the 
boundary condition near the secondary. For a fixed Tb c , 
Eq. (|28|l gives both F and D v , from which the surface den- 
sity and central temperature follow from Eqs. (|21H22|> . 



T c = 



8 2 / 5 K/fc) 4 ^ 1 / 5 pv-w* 

( 97r )3/5 a 4/5 K l/5 ^4/5 
(8tt) 2 / 5 Ql/B 



Km 2 + T hc ^ 



3/5 



(29) 



2 n 

A/fi 2 +T bc — 



-/"> 



(30) 



where 

^(0+4)710 (fi % f/^ c [k/^m p )] 2 / 5 ft 9 ' 1 



(31) 

Thus, solving the disk structure in these zones amounts to 
finding the torque at the boundary, Tb c . 

If T v {r m in) = then Eq. (JTTJ) shows that, 

T hc = -Mr 2 min n{r min )+ [ d r T d dr . (32) 



In practice, r m i n = risco for a disk without a cavity, and it 
is the inner edge of the disk if it has a cavity. Depending 
on which term dominates in Eq. (|27p . we distinguish the far 
zone (|T bc | < Mr 2 Q.) and the middle zone (|T bc | > Mr 2 Q). 
The far zone can be either well inside or far outside the 
secondary's orbit, but the middle zone is always outside. 
Well inside the secondary, the second term can be neglected 
in Eq. (|32p . and well outside of it, the second term dominates 
and the integration domain can be extended to oo. In both 
cases, The is independent of r. 

Eqs. (|27p an d (|32p show, that in the region outside the 
secondary, Tb c represents a torque barrier due to the sec- 
ondary's tidal effects. This parameter can also be used to 
obtain the migration rate of the secondary. Indeed, combin- 
ing Eqs. © and gives 

2Tbc 2Tgw 



m s r s tt s m s r a Q s 



(33) 



4-1.1 Far zone - unperturbed disk 

Without the secondary Eq. (|12p implies that Tb c = 
-Mr 2 sco ^isco. Substituting into Eqs. (f27H28)l , gives D v 
and F. P lugging int Eqs. flU), ([To), a nd (|21ll22p leads to the 
standard IShakura fc Sunvaevl (|l973h solution 



So 
F 

VrO 

H 



, 7 „ , ft 5 g -4/5 . 3/5 ,,1/5 

4.7 x 10 ~ce_-, m -i M 7 r 

cm 2 



3/5 3/5 

2 "r 



-VB • 2/5 .,-1/5 -9/10 2/5 



5.4 x W°KaZi< mJ 1 J M 7 ' r 2 ' ip 
7.9 x 10 13 -^m_iMf 1 



—MO 2 
8tt 



3 

r 2 f 



(34) 
(35) 
(36) 
(37) 



1.5M. m_i(y3 if P < 1 , 

0.28M. Q- 1 / 10 mL / 1 5 M 7 - 1/10 r 21/2 V 1/5 if ~ 1 . 

(38) 

Here and below, the subscript o denotes quantities related to 
the unperturbed disk, a~i = a/0.1, rh-i = (Af/AfEdd)/0.1, 
M^dd is the Eddington accretion rate for 10% radiative ef- 
ficiency, g_ 3 = q/10 -3 , M 7 = M./W 7 M , r s2 = r s /10 2 M., 
and we introduced 

<p=l — r I 2 sco fiisco/(r 2 fi) = 1 - (risco/r) 1/2 . (39) 

Without the secondary, in the radiation pressure dominated 
regime (j3 <C 1) the scale-height is approximately constant, 
and increases approximately linearly further out where gas 
pressure dominates (p ~ 1). 

The viscous torque, for future reference: 



Mr 2 9.if = 7.1 x 10 47 ergmM 2 r2 /2 v?, 



(40) 



where r 2 = r/10 2 M.. Sufficiently far from the secondary, the 
disk is independent of the secondary and follows Eqs. (|34F 
I38p with tp ss 1. However, the disk structure depends on the 
rate at which gas is allowed to flow in through M. 

4.1.2 Middle zone 

Now let us consider the opposite limit, T^ c S> Mr 2 Q, where 
the steady-state perturbation to the torque is significant. In 
terms of the dimensionless torque barrier, 



k 



The 
Mr 2 Q, 



r 2 A 
' r 2 Q, 



(41) 
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the formulae describing the unperturbed disk, Eq. (|34H38[) . 
get modified by replacing the boundary term with 



tp — > 1 + k . 
The disk quantities change to 



3/5 



Eq oc r 



-Tin 



2/5 



T c0 oc r 



-9/10 



-11/10 



H m = 



k + 1 

¥> 
fc + 1 

fc + 1 

¥> 
fe + 1 

(fc + l)/f /¥> oc r~ 1/2 



F oc r" 7/2 

-3/5 

V r Q OC r 



-1/10 



if P < 1 ; 



((fc+l)/co) 1/5 /f ocr 19/20 ifj9~l. 
and the migration rate follows from Eq. (|33|l 



2Tb c 



-2k 



Mr s 



(42) 

(43) 
(44) 
(45) 
(46) 
(47) 

(48) 



m 8 r 8 O s m s 
where we have assumed Tgw -C Tbc- Here and below, the 
superscript "m" labels the middle zone. Note that the di- 
mensionless angular momentum flux k can be interpreted as 
a brightening factor in the middle zone relative to the unper- 
turbed disk; k a is representative of the maximum brighten- 
ing, if k(r) is extrapolated to r s . In practice, the maximum 
brightening is even larger than k s in the near zone due to 
tidal heating (see § 14.31 below). 

The disk is modified within a radial range where the 
dimensionless angular momentum flux satisfies k > 1. This 
sets the outer boundary r™ of the middle zone, where the 
disk transitions to the far zone. From Eq. (|41[) . 



GM.M 2 



(49) 



Eqs. (|43H47[) represent a disk with negligible inflow of 
angular momentum but an inner boundary condition with 
a large viscous torque, corresponding to the torque barrier. 
Such solutions are oft en (somewhat misleadingly) refe rred 
to as a decretion disk (|Pringlelll99ll ; lLodato et al.ll2009h . To 
avoid confusion, we emphasize that there is accretion (i.e. in- 
flow) in this region, too, with a fixed M. However, the radial 
accretion velocity is greatly reduced, while the surface den- 
sity, temperature, and scale-height are all greatly increased, 
relative to an accretion disk around a single compact object. 

So far in this subsection, we have derived a solution for 
an arbitrary torque barrier or k, without specifying its value. 
In general, k is given by Eq. (|32p . which depends on the tidal 
torque in the near zone. Thus, to complete the derivation of 
the disk structure in the middle zone, we are first required 
to obtain the disk structure in the near zone (which we will 
do in § 14.21 below) . However, in the case of the steady-state 
cavity, the particular form of the boundary condition allows 
us to directly infer k, independently of the near zone, up to 
a factor A of order unity, which we show next. 



4-1.3 Middle zone - steady-state disk with a cavity 

When the tidal torque is sufficiently strong to clear a gap so 
that the secondary and the nearby gas move with a similar 



velocity, T v ~ The can be substantial over a large range of 
radii. From Eq. (|13H14|I . this requires 



rpUlg _ mlllg _ 



Here and below, the superscript 9 refers to solutions with a 
gap, and r g = Ar s is the outer radius of the gap. For this 
value of Tbc, Eq. (|28|l gives F and D v , and E follows from 
(f2l"jl . However, since the right hand side (RHS) of Eq. ([50)) 
depends on E itself, this gives an algebraic equation for Tbc- 
The solution is 

i 3 / 4 ,. 3 / 4 

(51) 



Q3/4 1/2,1/8 



47r 1 / 4 ( M m p /fc) 1 / 2 cr 1 / 8 
= 1.6 x lO^erga^m^A- 11 / 16 ^ 8 ^/^ 4 ^ 3 / 8 . 

Note that this is independent of the tidal torque model (i.e. 
the A or / in Eq. ©), since here the tidal torque is set by the 
boundary condition of the gap. This solution breaks down, 
and becomes tidal torque dependent, if the gap closes, which 
we discuss in § 14.21 below. 

The dimensionless angular momentum flux from 
Eq. gTJ is 

fc ms = 23ai / 1 2 m ;? /8 Af- 3/4 q% 8 A- 19 / 16 r- s2 7/8 ( 

(52) 

In particular, near the secondary fc mg (r- s ) = m s /(47rr|E g ) 
is the ratio of secondary mass to the accumulated local gas 
mass The only free parameter in this zone is A, which we 
determine explicitly in § l4.3.2l below. 

In the range r s < r < r f mg , fc mg > 1 and Eqs. (|43H47)1 

give 



E mg = 3.1 x 10' 



• 3/8 3/8 
m -l 1-3 ^-33/80„-9/40„-9/10 



cm 2 Q V 2 M 1/4 



1/4 



T" g = 1.9 x itfKmlt-^TE^ 11 ' 



40 -3/20 -11/10 



-1/10 



(53) 
(54) 

(55) 



{ 35M. aL / 1 2 m 5 _ / 1 8 M 7 - 3/4 A' 19 / 16 g% 8 r s2 7/8 
x(r/Ar s )-V 2 if /3 < 1 , 

0.53M.mi /8 M- 1/4 A- n / 16 ^ / 3 8 r s - 3/40 rf /2 
if p ~ 1 . 



(56) 

where E(r) ~ at r Ar s . Outside r > r" lg , fc mg w 0, 
and the disk approaches the unperturbed solution given by 
Eqs. (|34H38[) . In the transition zone, between the middle 
and far zones, r ~ r™ g , one needs to use Eqs. (|43H47jl with 
k = fc mg (Eq.[52j. 

When a cavity is present, the migration speed of the 
secondary follows from Eqs. (|33p and (|51|l : 



-550- 



Va^S/S A/f l/4^-3/8 A -ll/16^1/ 



(57) 



14 Here fc mg (r s ) = B 5 / 8 using the Syer-Clarke parameter B = 
m s /(47rr 2 E s0 ). 
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This expression is consisten t with the secondary d ominated 
Type-II migration rate of ISver fc Clarke! (|l995t ) who as- 
sumed A = 1. Note that the migration speed is slower 
than the gas inflow velocity without the secondary, |n^ g | < 
|iVo(t"s)| m Eq. (J37J). This is referred to as disk-dominated 
Type-II migration, which is appropriate if the secondary 
mass is smaller than the unperturbed local disk mass m s ^ 
47rrgEo(r g ) (or equivalently fc mg ^ 1), but large enough to 
open a gap. 

It is interesting to note that the structure of the middle 
zone does not depend explicitly on the tidal torque model, 
d r Td (in particular A or the fa parameter in Eq. [5]); the de- 
pendence is implicit and arises only by fixing the value of 
A. Physically, while the tidal torques are negligible in this 
region, the effects of the tidal torques are still communi- 
cated to the region by setting an effective hydrodynamical 
boundary condition. We determine A in § 14.3.21 below and 
find that, in fact, it only weakly depends on d r T^. 



4-1-4 Consistency of steady-state 

A basic assumption of our model is that the radial struc- 
ture of the disk is in a quasi steady-state as the secondary 
migrates slowly inwards. To check the consistency of these 
steady-state solutions, we must verify that the implicit time- 
dependence in the surface density profile through r s (t) does 
not violate the continuity equation ([TJ significantly, so that 



the overflowing solution over a wide radial range. For global- 
steady state, a necessary condition is 



d r M = -d r (2TvrT,v r ) = -2nrd t T, = 



(58) 



Integrating over radius, the relative error in the accretion 
rate 



/27rrE mg dr _ / 27rri; s r ? g d r . s £ mg dr _ 9 



M 



27mv lg E m s 44 V Xr, 



11/10 



(59) 



which is ~ 20% near the gap edge. The error in Tb c based 
on Eq. (|51|l is ~ 13%. However, the error in the accretion 
rate exceeds unity at large radii, outside 4.2Ar s . The steady- 
state assumption breaks down because as the secondary mi- 
grates inward, the steady-state gas density near the edge 
of the gap continuously increases with time. If M is fixed 
near the gap edge to be a constant fraction of the Edding- 
ton value, the true accretion rate M(r) at larger radii must 
be larger, to supply material for the increasing gas density. 
Conversely, if M(r) is fixed at large radii, then it becomes 
smaller approaching the gap edge. Such no n-st eady-state so- 
lution s have been derived bv lPrinele] (|l99ll ) and llvanov et al.1 
|l999l ) by solving the nonlinear diffusion equation ([7| for a 
fixed outer boundary condition, assuming that the viscos- 
ity can be expressed as v = fcEV 6 where fc, a, and b are 
constants. In particular, llvanov et al. I i|l999h derived a non- 
steady, but self-similar solution. In that solution, the migra- 
tion is slower, an d the angular mome ntum flux is lower, 
compared to the ISver fc Clarke! (|l995l ) steady-state solu- 
tions with a fixed M f or the same binary and disk parameters 
(|Sver fc Clarkdl 19951 ). The quasi-steady migration rate and 
brightening factors f or a truncated disk wi th A > 1 are inter- 
media te between the lSver fc Clarke] (|l995l ) and llvanov et al] 
(|l999h solutions. 

The steady-state condition is typically not violated in 



r o r f° 27rrE m dr 




7s s 


M 




2 + 7Sr 



fc. 



19/ 



47rr s 2 E (r s ) 



< 1 (60) 



where we have used Eqs. (|29|) . (|48|) . and (|49|) and defined 
7e s = 91nE m /ainr s and 7& = 91nE m /91nr. This sets 
a maximum limit for the dimensionless angular momentum 
flux fc s . For a /3-disk, this implies 

fcsmax = 3.2 ^I^V^ml^Vr^ 19 ^ 7719 , (61) 

and fc sm ax is larger (i.e. less restrictive) for radiation pressure 
dominated ev-disks. 

An important qualitative differe nce between th e over - 
flowing model presented here and the lSver fc Clarke] (|l995l ) 
model for a truncated disk is that 7s s > for the former 
as we show below. In contrast, in the overflowing case, the 
excess surface density and the dimensionless angular mo- 
mentum flux fc a in the middle zone both gradually decrease 
during the inward migration of the secondary. Thus, the ex- 
cess surface density diffuses radially outwards. If fc s < fcsmax 
then the diffusion is sufficiently fast to reach a global quasi- 
steady-state throughout the middle zone. If this is not satisi- 
fied, then the outer parts of the middle zone cannot respond 
as quickly as the object moves inwards and the structure of 
the disk in these regions will depend on its previous history. 
However, since the viscous timescale is always much smaller 
than the migration timescale in at least the inner parts of 
the middle zone, the local disk structure of the overflowing 
solution might approach an approximate steady-state there 
with a constant M even if fc s > fcsmax- The migration rate of 
the secondary depends on the near-zone of the disk, which is 
expected to remain insensitive to perturbations in the outer 
parts of the middle zone in an overflowing disk(3 We leave a 
detailed investigation o f the time depend ent overflowing so- 
lutions to future work (|Salem et alj|20l"3 . in preparation). 



4.2 Near zone 

Now let us consider the regions near the secondary where the 
tidal torque and heating are important. We discuss steady- 
state solutions inside and outside of the secondary's orbit, 
in turn, without and then with a circular cavity. Deriving 
the physical properties of the disk in this region is useful 
to provide an estimate of the torque barrier, The, at the 
interface between the near zone and the middle zone. As 
explained previously, the torque barrier sets the overall scale 
of the physical properties in the middle zone, as well as the 
migration rate of the secondary. We therefore first compute 
the value of the torque barrier for an overflowing disk, as 
well as for a disk with a wide gap. In the latter case, we 
then compare the value with the torque barrier in the middle 
zone derived above (Eq. I51|l . By equating the two, we can 
estimate the gap size (i.e. A), and obtain the conditions for 
gap opening and closing. 



15 This is different from a transient truncated circumbinary disk 
where the migration velocity is comparable to the local gas accre- 
tion velocity, which can exhibit hysteresis throughout the middle 
zone llRafikovll20l3) . 
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4-2.1 Inside the secondary orbit 

Consider the region just downstream the secondary, outside 
the torque cutoff A = r s — r > H in Eq. (|6]), assuming a 
steady-state overflow (i.e. no hollow circular cavity). Based 
on Eq. the viscous torque decreases in the vicinity of 

the secondary and for a sufficiently large secondary mass, 
the angular momentum exchange is dominated by the tidal 
torque. In this regime, 



\d r T„\ « \d r T d \ = 27rr|A|E w Md r (r 2 n) 



implying that 



4tt |A| 



M 



1 A 



2tt fq 2 r 2 Q 



(62) 



(63) 



for a Keplerian disk (here and below, the superscript m refers 
to the solutions in the inner near zone). Eq. (|20l) shows that 



F ni « D d = i(fi s - fi)AE = - fi.) -> § - . (64) 

2 87T 2 r 



The asymptotic limit corresponds to A < r. From Eqs. ((8} 
and 115} 



2A , 2 r 5 Q 
rti = fq "a*" 



km n - n s 



#o A 

-)■ — if j9 ■ 

2 r 



(65) 
(66) 



The latter equation shows that the secondary makes the disk 
thinner downstream if the disk is radiation pressure domi- 
nated. Combining Eqs. (|63H64[1 with (J2TJ, gives D™. The 
viscous torque then follows from Eq. (|18[) . To first beyond 
leading order, for 6=1, 



E nl = 



: 5.7 x lO^/I^m-i^ra 1/a (A/r) 4 , 
T" 1 = 1.5 x 10 6 K fZl ,A rnl?Mi 



-1/4 -1/2 -7/8 
2 



5/16 



5/4 



i^ ni = 3.9 x 10 13 J^-rh_iMf 3 f- 



s cm 

on Cm 2 -1/2/ A / 

= 30— f-2Q- 3 r 2 (A/r 



2 
-4 



5/4 A 
r 



f 0.75M.rh_i(r/r s ) 5/4 (A/r) if /3 < 1 . 

# ni = 0.36M./I a 1/B m^ 1 4 M 7 - 1/8 9 Z 3 1/ VJ r/16 
i x(r/r s ) 5/32 (A/r) 5/8 if /3 ~ 1 . 



(70) 

(71) 

(72) 
(73) 

(74) 



During the inward migration of the secondary, r B decreases, 
and the surface density evolves in a self-similar way. The 
surface density, midplane temperature, surface brightness, 
and scale-height all decrease significantly near the secondary 
with large qs- The radial flow velocity becomes very large 
in the close vicinity, and the flow may become advection 
dominated there. However, we do not extrapolate this so- 
lution inside the Hill radius of the secondary because the 
torque model is invalid there. 

It is remarkable that for a fixed M, the fractional 
perturbation to the surface brightness and the radiation- 
pressure-dominated scale-height are universal in this region, 
independent of the binary and disk parameters. This prop- 
erty is general for an arbitrary torque or viscosity model in 
radiatively efficient steady-state disks. The surface density 
in this regime is also independent of the viscosity model but 
it is sensitive to the torque model: it is set to ensure that 
the tidal torque matches the angular momentum flow asso- 
ciated with M. The original value of the surface density is 
suppressed by a factor proportional to g~ 2 A. These solu- 
tions are valid only for disks with relatively large secondary 
masses, such that A n i > j"h, but in which there is gas inflow 
across the secondary orbit. 



T-ini r, a . ln 50 ■ 3/2 ,.-5/4 „ ,7/4 -5/2 

T„ =2.4x10 ergQ_im_ 1 /_ 2 M 7 ' q_ 3 

/ \ 5 /!6 / A \ 21 / 4 

T" 1 exhibits a sharp cutoff near the secondary. 

We can verify that the working assumptions hold in this 
region. Eq. (|66[) shows that A > H holds for all A, since the 
unperturbed disk is thin, H a < r. Since oc T°* oc A 21/4 
which implies that D" 1 F nl is indeed satisfied for suffi- 
ciently small A. Coincidentally, the assumption in Eq. (|62|) 
is satisfied within a distance A n i from the secondary, where 

rs 1 + Y^x n [ 

and 

x ai = 0.1 aZt^rh-^M^f'J^q 1 ^ 17 ^ 34 . (69) 
The disk parameters in the region r s — A nl < r < r s — ru 



4-2.2 Outside the secondary - unsaturated torque 

Next consider the region just outside the secondary. Here we 
examine the case where the secondary is massive enough for 
the tidal torques to be important. After a significant amount 
of gas pile-up the viscous torque eventually counteracts the 
tidal torque and creates a stationary inflow. In this regime 
the tidal and viscous torques counteract one another and 
both greatly exceed the momentum flux in Eq. (|ll|l such 
that 

Md r (r 2 Q,) < d r T v w d r T d = 2-ivrAY, . (75) 

The tidal heating rate is much larger than the viscous heat- 
ing rate of a disk without a satellite, making the disk much 
hotter and thicker. Let us assume Dd 2> D v here, Eq. (|20|l 
implying that 

F^D d = i(n s -fi)AE. (76) 

In this subsection we examine the case where the torque is 
not saturated, H(r) < r — r s , so that A = r — r s in Eq. (|5J). 
This is most relevant for relatively small mass ratios, (i.e. 
typically q < 10~ 3 , sec Eq. l95l below and Paper II), where the 
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banking up of the stream in this region is modest. Substitute 
Eq. (JUJ) for E with b = 1, and use Eq. JTSJ, 



o \ 4 /5 



Q.)Ar 



-2/5 



5 J 1 



4/5 



7?l/5 



(77) 



where as is the constant coefficient in Eq. (|2f |) . Solve this 
for F and plug back into Eq. (|76[l 

a r T„ _ 2 13 / 19 tt 13 / 15 rA 5 / 6 

rj /3 ~ 3 2 /"4 /6 (n. - fi) 1 / 6 ( j 

Integrating both sides between n, the inner edge of this 
region (see further below for a discussion of the value of r\) 
and r, 



rj /s (r)-r^ 8 (n) = 



, A 5/6 



dr. (79) 



We are most interested in the case where the tidal torques 
increase T„ substantially in this region. If so, we approximate 
T v (n) = 0. Substituting Eq. ([5}, and rearranging gives 

2T U = 3.5 x lO^ergallM^V^^s^^C'^r^n) (80) 

where the superscript neu refers to the case of unsaturated 
torque in the external near zone, and 



c 



•As 



-1/6 



(81) 



n/fs 

2\ 1/6 2(n/r s )- 115 / 72 



5 (A;/r s ) 5 / 2 



1 - 



115/72 



5/2' 



In the second line, the approximation is accurate to within 
6% for Ai = n - r 8 < 0.3 r 8 . Note that T^ cu depends on 
radius only through (; it increases monotonically and ap- 
proaches a constant value, which depends very sensitively 
on Ai/r s . In practice, one might expect 



Ai ~ th and n = r s + th 



(82) 



for a disk without a cavity because the tidal torque model 
is valid only outside this region, and within this distance the 
gas may flow across the secondary orbit along radial streams 
or horse shoe orbits. In the following we keep Ai/rH general. 
We incorporate a factor of (rn/r s )~ 5 ^ 2 in the prefactor of 
Eq. (|80|l and introduce a renormalized £, as 



Cr 



5/2 



1/6 



C 

2(n/r s )- 115/72 



(83) 



3 J 5 (Ai/r H ) 5 / 2 
Then Eq. (JSOj) becomes 

„ = 1.7 x 10 erga_ l M 7 ' /_ 

115/72 



r; \j 115/72 

r 



5/2 



/2 9- / 3 2 r- s 1 2 /4 CR 

A^ 5/21 3 
A 



(84) 



where 



lim T™ u (r,r s ,n) 



^. \ - 15 / 2 / A. \ -H5/24 



x 



1+* 

r s 



(85) 



(86) 



The outer edge of this region is where Eq. (1 751) is first 
violated, i.e. where the viscous torque density! 16 ! becomes 
comparable to the accretion term d r T„ ~ Md r (r 2 Q). We 
substitute <9 r T„ from Eq. ([78)) utilizing (|80} and get 



1 = 



d r T" 



Md r {r 2 D.) 



no. in— 7 -2 --1 71 ,,1/2 ,5/2 5 -1/4 

= 2.8 x 10 a.jm.iM/ /_' 2 g- 3 r s2 



-11/24 



-7/2 



((r,r s ,n) 2 



(87) 



We label the radial distance of this interface from the sec- 
ondary as A£f u . While this equation of a single variable can 
be easily solved numerically for A£f u for any fixed r s and Ai, 
we may derive an analytical approximate solution as follows. 
Assuming 1 S> AJ}f u > 3Ai, £ is close to its asymptotic max- 
imum, which implies 



A, 



5.0a- 4 / 7 m: 2/7 A/ 7 1/7 /- / 2 14 g3 13/2 V s " 2 1/14 



[1 + (A ; /r a 



-115/126 



(A ; /r H ) 17 / 7 



(88) 



To obtain the dependence of the physical parameters 
on radius, we first derive D v and F by substituting Eq. (|84[) 
into JTSJ and |(77]). Then E and T c follow from Eqs. (J2TJI and 
In the range n < r < r s + A^ cu , 



1.0 x 10 7 - 



cm 

'A\ 1/2 
I -J Cr, 



-2 ,3/2 ,,1/2 4/3 -3/4 



-23/24 



^ c 



-25/24 



6.8 x 10 



12 erg 



-2 ,5/2 10/3,,- 

a_i/_ 2 g_ 3 M 7 



1/2 -13/4 



-77/24 



5/2 



<,R > 



,™ Cm 2 . ,-3/2 ,,-1/2 -4/3 -1/4 

173 — a_im_i/_ 2 M 7 ' 



g_ 3 ' r s2 — 



-1/2 



if™ = ^ 



Cr 2 , 

2 ,5/2 Ml/2 -1/4, 



(89) 

-1/2 

Cr- 
(90) 

(91) 

-1/24 

(92) 



f 0.13M. aZlfJiM^r^'^r/r.)-^ 

x(A/r s )- 5/2 Cl ifj9«0, 

0.23M.a:; /2 /y 2 2 r s2 (r-/r s ) 47 / 48 (A/r s )^/ 4 C R /2 
if /3 « 1 . 

(93) 

Here r = r s + A, and these equations are formally correct to 
first beyond leading order in A/r a , but we find them to be a 
good approximation typically within 15% even for A/r s > 1. 
Interestingly, all of these physical parameters have a local 
extremum in this zone. We label the distance corresponding 
to the maximum local disk luminosity 47rr 2 F(r) as A°™ k . 
We find that Ap°" k /Ai is a slowly decreasing function of Ai, 
which varies between 1.55 and 1.4 for < Ai < r s . 



16 Note that matching the derivatives at the interface does not 
contradict T£ eu > Mr 2 !i there. 
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To match T„ at a radius r = r s + A[^ u . the interface 
between this region and the middle zone, we must set Tb c = 
rr u (rr). Based on Eqs. (J84J) we may use 

if AST » A ; . 



rj-iU 



(94) 

This equation is typically valid in the overflowing case (see 
Eq. I88|) . as long as q is large enough that T"m ax > T^o^s) 
(strongly perturbed solution), but not too large so that H < 
r — r s . We discuss the solutions if the latter condition is 
violated in § 14.2.31 below. Comparing Eqs. (|40|) and (f86)l 
shows that the minimum mass ratio to cause a significant 
gas buildup with unsaturated torques: 



ncu 
*7mm 



9 x 10 



-4 4/5 • 2/5 



,V7- 



1/5 ( -l 1/10 



2 ' S 2 



(95) 



For smaller masses, the disk structure is not modified signif- 
icantly outside the secondary, and one has to use Eq. ()67|) 
for the torque at the inner boundary in Eq. ()79|) with 
T v [rC) = T"'(r — ru)- We do not show these more general 
but more complicated expressions here. 

We note that the results in this section are sensitive to 
Ai if different from ru, which sets the distance at which the 
gas can flow in across the secondary's orbit without signif- 
icant resistance. While Ai ~ ru is reasonable based on the 
horse shoe orbits in the restricted three body problem, we 
keep it as a free parameter in the following. 

4-2.3 Outside the secondary - saturated torque 

Here we again assume that Eqs. (J75H7BJ hold, but now ex- 
amine the case of much higher secondary masses, where the 
scale-height is increased so much that H(r) > r — r s and the 
tidal perturbation enters the torque cutoff regime. We make 
the simplifying assumption here that H(r) > r — r B holds 
throughout the near zone so that the tidal torque in Eq. ((5} 
does not alternate between saturated and unsaturated. It 
is straightforward to obtain more general solutions, but we 
find this exclusively saturated OR unsaturated assumption 
to be an excellent approximation in most cases. 

Due to the large torque barrier and tidal heating, the 
disk in this region is typically radiation pressure dominated, 
and we accordingly assume p ~ for the analytical solutions 
below. 

Eqs. JT6J, GM, and (ED show that H = ^-Q,~ 2 F and 



£ = asr 



8/5,^4/5^,-1/5 



for 6=1, where as is a constant. 



Substituting into Eq. (|76l) . 



(96) 



This equation can be solved for F as a function of r and T v . 
Plugging back into Eq. (|75[) leads to a separable first order 
differential equation for T v 



d r T v 

1 v 



air. 



9/52 



n s - n 



-21/26 



25/13 / f \ 14/13 



(97) 

where a\ is a constant independent of r and r s . Now 
usj^l fi/fl s ~ (r/r s )~ 3//2 and integrate both sides assum- 
ing T u (ri) w for some n > r s . Here n is the radius where 

17 The disk rotates with nearly the local Keplerian angular 
velocity, but slightly slower due to a radial pressure gra dient: 
(n K - n)/ft K oc H 2 /r 2 see Eq. (78) in lKocsis et""aH feulll) . 



the torque model breaks down, for which we adopt the Hill 
radius around the secondary n = [1 + (q/3) 1/,3 ]r s if a gap 
does not form. Thus, 



(98) 



T^\r)-Tl^\ ri ) = ii ai r s 61 / 5 V(r,r s ,n), 



where we introduced a dimensionless function 



. 2 „{Sk-(l 5 7\ 2 / O-jlg 5 7, 
^=3 B (^^ ; 2^13 ~3 B -^ ! 26'13 ( " } 



3.7 
0.4 In 



- r, \ 5 / 26 



r — r 
r 

r — r B 

r — n 



5/26 



(100) 
(101) 



Here B(x; a, 6) = Jjf t a ~ (1 — t) ~ dt is the incomplete beta 
function, and the last two lines are simple approximations, 
typically accurate to within 15%. We can now use Eqs. ()18|) 
to get D v . For T v (n) f» 0, after substituting the value of a\, 
Eq. (|98j yields 



rr = LI x lO^ergal^^/^g-^'Mr 5722 ^ 1744 ^ 13711 

(102) 



DT = 1.2 x lO 15 ^^ 11 /^ 22 ^ 1 ^- 21 / 22 
cm 2 s 



-93/44 / r 



-7/2 



,13/11 



(103) 



Here the superscript nes refers to the case of saturated torque 
in the external near zone. Note that T„ es depends on radius 
only through ip. Close to the inner boundary of this region 
n, it grows quickly with 8r/r and saturates to a constant 
at Sr/r ~ 1. This can be understood, since this solution 
neglects angular momentum flow, the viscous torque is equal 
to the integrated tidal torque density, and the latter has 
a cutoff at Sr/r ~ 1. We can verify that T,' lcs > T v0 is 
indeed satisfied in this region (c.f. Eq. I40p and so the first 
assumption, Eq. ([75}, and T v {ri) ~ are well justified. 

Can we use the asymptotic maximum of T" cs as an es- 
timate of the torque at the outer boundary of this region, 
to estimate T™ in the middle zone of an overflowing- disk? 
In many cases no, because the disk transitions to the mid- 
dle zone much closer Sr ^ r s , implying that the torque at 
the outer boundary of this region can be much less than its 
asymptotic maximum. The outer boundary of this region is 
where H — r — r s . To figure out exactly where this happens, 
we proceed to determine the disk structure in this region. 

Now Eqs. ffT03T> and J96} give D v and F; all other disk 
parameters then follow from Eq. ©, JT6J, and (|21H22|) . The 
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result within n ^ r < r B + H nes is 
S ncs = 2.6 x 10 6 ^alf /11 M 7 5/22 /! / 2 22 ^ / 3 11 r s - 2 3/44 



a - n 
a 



-1/26 



1.6 x 10 D Ka:f /11 Af 7 _z/11 /-2 

-25/26 



-27/26 



10/11 



(104) 



2/H ^l/H„2/ll r -6/ll 



a 

9.8 x 10 1 



3 ' s2 
3/11 



(105) 



■ ct_ 



a 



5/26 



lvl 7 7-2 9-3 T s2 

-73/26 



2/11 



(106) 



66 o£ n i Q «/ 11 m„ 1 M 7 - 5/22 /r2 3/22 g:3 3/11 r s - 2 41/44 



a - n 



1/26 



, 1/26 
r ^ V- 10A1 



19 M. al 2 / 1 ^/^ 22 /!/ 22 ^ 11 ^ 44 



a - n 
a 



5/26 



5/26 







2/11 



(107) 



(108) 



We can now confirm the consistency of the second as- 
sumption, Eq. ([T6]>. using Eqs. (|103j) and (|106[) . Indeed, 
^" cs < F ncs near the secondary since D° es scales with a 
higher power of "0. £ and T c have a maximum, iv decreases 
and becomes practically constant, while H slowly increases 
in this regime. 

The outer boundary of this region, r™ a , is wherj^l 



r rues / lies \ r n 

H (r m ) = Sr n 



(109) 



We use s as an approximation to the transition radius to 
the middle zone. After substituting Eq. (|108l) . Eq. (|109l) is 
a nonlinear algebraic equation for rjjf s . While it is easy to 
solve it numerically for any choice of parameters, it is still 
useful to derive approximate analytical solutions. We find 
the following method yields results that are accurate within 
20% for a wide range of parameters. Use Eq. (|100|) . expand 
H ncs to second order in Sr^ B , use (1 + ax) fts (1 + x) a for 
small x and 1 — x a ~ — alnx- for small a. This gives an 
approximate relation 

^S- » 0.17a:t /17 M 1/17 /! / 2 17 ^° 3 /17 r S 2 5/34 



f; Ar nos 
x [ A l n ^HL_ 

. 26 5n 



4/17 



(110) 



where Sri = n — r s . Eq. (|110p can be solved analyt ically 
using the Lambert W- function (jCorless et al.lll996lrn 



W-i(-a) « ln(a) - ln(-ln(a)) 



(111) 



Sri 



5r\ 



W-i(-a) 



4/17 



Sri exp 



17 



W-i(-a) 



(112) 



18 Note that in this section we use <5r = r — r s instead of A since 
in this region A = max(<5r, H) = H (see Eq. [BJ. 

19 The Lambert W-function is denned to be the inverse of the 
function /(W) = Wexp(W), where we need the real branch with 
the larger absolute value, W_ i, defined on / > — e _1 = —0.368. 
The approximation in Eq. J 11 1 D is correct to within 20% for all 
< a < 1/e. 



where the two forms are equivalent, and we have introduced 

n , C r ,,-1/4,-5/4 -13/12 5/8 {Sri \ ^ ^ 

a = 0.465 a_iM 7 f_ 2 > q_ 3 r s2 ' 1 — 1 
Finally, we substitute in Eq. (|102[) and use Eq. (|101[1 . 



(113) 



T™ s (r™ B ) = 6.9 x 10 47 erga:^ /1I A/ 7 ' o/ ^7_ 2 



2/H (1^45/22 f 5/22 5/11 61/44 
y-3 's2 



x [-W-i(-a)] 13/n 



(114) 



Note that W-i(— a) depends logarithm ically weakly on 
the disk parameters (jCorless et al.l fl996l ) . in practice 1 ^ 
|W-i(-a)| < 10 for 10" 4 < a < 1/e = 0.368. Here a < 1/e 
is required for this solution to exist, implying that q and r s 
have to be sufficiently large and small, respectively. In the 
opposite case, the torque is unsaturated (§ 14.2.21) . 



4.3 Transition between near and middle zones 

4-3.1 The case with overflow 

The value of T v at the outer edge of the near zone, is to be 
matched with that in the middle zone, Tb c . If the tidal torque 
is unsaturated in the near zone, we approximate Tb c with 
the asymptotic maximum value, T^ c = T°™ ax . Otherwise, if 
it is saturated, then we set T bc = T" es (r£ 1 es ). Matching the 
middle and near zones at r^ s assumes that T„ does not grow 
substantially in the transition region between the saturated 
near zone and the middle zone, i.e. outside of rjjf 8 but within 
a radius where the tidal effects are still non-negligible. We 
find this approximation to be better than 10%. Thus, to 
match the torque at the outer boundary of the near zone 
and the inner boundary of the overflowing middle zone, we 
combine the saturated and unsaturated cases as 

min{T- s (C cs ), T,%i u ax} if a < 0.368 , 
T™, if a ^ 0.368 , 



(115) 



where T? ea (r™ a ) and T^ ax are given by Eqs. (ggj and {TT4J). 
If this satisfies T b ™° < T™ 8 for Sn = r H (see Eq. [STJ, then 
the satellite migration velocity is less than the gas bulk local 
inflow velocity, and the overflowing steady-state solution is 
self-consistent. In the opposite case the disk forms a gap 
with Sri > ?"H. 

The dimensionless angular momentum flux in the mid- 
dle zone (Eq. QT} is 



A" 



1.3 a: 



-1/4 



1 + 



(I) 



1/3 



-15/2 



0.97 a 



-2/11 



m_ 1 M 7 < f. 



1/22 ,5/22 5/11 39/44 



|W(a)| 



(116) 

13/11 

(117) 

Gap overflow causes the torque level to decrease in the mid- 
dle zone, which suppresses k. We discuss gap closing in § 14.41 
below. 

The disk flux in the near zone is larger, due to tidal heat- 
ing, than in the middle zone. To show this, we next compare 
the luminosity of the near zone to the middle zone explic- 
itly. In the case of unsaturated torques in the near zone, 
we find that the local disk luminosity, L ncu (r) = 4 7 rr 2 J incu 
(see discussion following Eq . 191 [I . peaks sharply near »"peak = 
r s + 1.5 th- We find that the integrated flux from within a 
ring of width 0.5 ru is approximately 



(118) 
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In the unsaturated case, rn/r a = (g/3) 1//3 <C 1, and so the 
net luminosity of the middle zone exceeds that of the near 
zone. For saturated torques, the maximum brightness cor- 
responds to the outer boundary of the near zone, r^ cs . We 
assume an effective radial width 5rJJf s given by Eq. (|110p . 
From Eqs. (|16[1 and (|109|l . assuming a radiation pressure 
dominated near zone, the dimensionless angular momen- 
tum flux of the near zone relative to the unperturbed disk 
is 

, ncs _ iirr^Sr^F^ir™) _ SrT H ncs (rT) _ (5r™ s f 



47r(rgf s ) 2 F (ri) 



H 



0.45 a 



-8/17„ 



(119) 

,8/17 



■ 1 /- J ' 1 ''JWr''^ r V,7"[-W(a)] 

Comparing Eqs. (|117|l and (|1 19[) . and recalling that typically 
1 ^ |W(— a) | < 10, we conclude that the luminosity of the 
saturated near zone can exceed that of the middle zone by 
a factor between ~ 3 — 10 for q ~ 0.1. 

Given the dimensionless angular momentum flux in the 
middle zone, the disk parameters and the migration rate 
are given by Eqs. (|43H48|) . We substitute Eqs. (|116H117[) 
to obtain an explicit formula for the disk parameters. The 
migration rate in the overflowing disk with saturated and 
unsaturated torques is, respectively, 

-115/24 

f 5/2 ~ 3/2 M 3/2 r 3/4 



cm 3 o 

— 30 ct_ 1 j_ 

s 



1 + '§ 



1/3 



cm _ 
23 — a 

s 



2/11 
1 



j.5/22 
f-2 1. 



-6/11 
-3 



,23/22 83/44 r 



Al 



(120) 
-W(a)] 13/n . 

(121) 



4-3.2 Disk with a cavity 



Let us next turn to the case with a gap. We determine the 
radial distance to the outer edge of the gap, Ar s , here by re- 
quiring that the migration velocity matches the rescaled gas 
inflow velocity in the middle zone as stated in the boundary 
condition, Eq. fTS) . With A in hand, Eq. (|53H57p determine 
the disk parameters in the middle zone. The solution is dif- 
ferent when the tidal torque is unsaturated near the inner 
edge and when it is saturated, which we discuss in turn be- 
low. 

First, assume that the torque is unsaturated all the way 
outside of the gap (H ^ r — r s ). We use Eq. (|87p to obtain 
C(Ar s , r s , n) at the interface between the near and the middle 
zone. We substitute in (|92[) to obtain the gas velocity at the 
interface rescaled by 1/A: 



tf ou (Ar„) 
A 



on Cm -f 2 -1/2-, -3/2,<-> ,\- 
30— /-2<7- 3 r s 3 'A (A - 1) 



(122) 



To obtain the migration velocity, we identify Tb c = 
r" cu (Ar s ) in Eq. and substitute Eq. (J80j), and elimi- 

nate £ using Eq. ((87)) . This gives 

v Br = 8022^ a_ 1 m 3/ 1 2 M 3/4 /r 2 5/4 g: 3 7/2 r s 9 2 /8 A 11 / 16 (A-l) 21 / 4 . 

(123) 

The boundary condition, Eq. (|13[) , states that Eqs. (|122[) and 
(|123[) must be equal. This provides a nonlinear equation for 
A. We solve this equation perturbatively. To first beyond 
leading order, 



A u = 1 + 5 U (1 + 5 U ) 



-35/148 



(124) 



where 



4/37 ■ -6/37 J-9/37 22/37 , ,-3/37 -13/74 



S u = 0.55 cCr 'm- 0/ *'fl' 2 ° qlT M 7 



(125) 



The 'u' subscript is introduced to distinguish the case with 
unsaturated torques. 

Next, consider the case of the torque cutoff. The for- 
mulas in § 14.2.31 are not limited to the overflowing case, as 
long as the torque is saturated (Sr = r — r E ^ H). If a gap 
opens then Sn marks the distance to the edge of the disk 
in Eq. (|98H99|) . f° r which 5r-[ > ru- Here, Sri can be elimi- 
nated using the boundary condition Eq. (|13p as follows. We 
identify A = r™ s /r s in Eq. (JT3J where r^ s marks the edge 
of the near zone according to Eq. f|109f) , so that A = 1 + 5i 
where <Si = 6r™ s /r B . Combine Eqs. (fT07tiT09]) to eliminate 
tp(ri) from the bulk gas velocity at ri 



3 (C CS ) nK cm. 2 -3/2 [!-(! + &)■ 
A = 0.15— m_i/-2g- 3 r s2 



-3/21 



(126) 
eliminate 



2.0 x W^a-.rJU-T^^ 5 ^ 5 ^ 



Similarly, assuming T bc = T™ B (r™ B ) in Eq. 
V>(r£T s ) from Eqs. (fT02]) and (TTMJ, we get 



82 [l-(l+(Ji)- 3 / 2 ] 5 /4- 

(127) 



The boundary condition, Eq. (|13p . states that Eqs. (|126 
and (|127l) are equal. After rearranging, we get 



r r\io -4/37 • 4/37 ,.9/37 22/37 „ ,-3/37 -33/74 

5i = 0.13 (*_/ m.'i f_' 2 q_i M 7 r s2 



3 + 3<5i + Sf 



2 T9/37 



l+(l+<$l) 3 / 2 



(l+<5l)- 17/74 . 



(128) 



The last two terms can be omitted within 10% accuracy for 
< 5 < 3. This gives the characteristic gap scale in the 
torque cutoff zone 



1+0.13 aZt ,37 m_!l 7 M^i°' f°J£ ' q"i a ' r& 



3/37r9/37 22/37 -33/74 



(129) 



provided that A s — 1 > rn/r s (i.e. r — r a > r s + fH); otherwise 
no gap is possible. 

Thus, the dimensionless angular momentum flux fol- 
lows after substituting into Eq. (|52p . In the unsaturated and 
saturated cases, 



fc s mgu = 23 a 
fc" lgs = 23 a 1 



1/2 . -3/8 
-1 m 0.1 



m: 



3/4 5/8,-11/16 -7/8 



l/2^-3/8 M -3/4^5/8 ^-11/16^-7/8 



1 '"-o.i 



a; 



(130) 
(131) 



The migration speed of the secondary in case of a gap with 
unsaturated and saturated tidal torques, respectively, is 



wf r u = -550— a'm^M'^q 
s 



1 / 2 ^,5/8 ^1/4^-3/8^-11/16^1/8 



1 / 2 ^,5/8 ^1/4^-3/8^-11/16^1/8 



= -550— a^th^M^q. 



(132) 
(133) 



These estimates depend on the somewhat arbitrary def- 
inition of A that we have adopted in the two cases. In the 
unsaturated case, we have identified it with the outer edge 
of the transition region between the near and middle zones, 
where d r T v — Md r {r 2 Q,), while in the saturated case, we 
considered it to be the inner edge of the transition region, 
where H — r — r s . While these conventions could be mod- 
ified, they do not affect the overflowing solution. They do, 
however, influence the secondary orbital radius where the 
gap closes, which we discuss next. 
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4.4 Gap opening and closing 

The previous sections define the disk uniquely, which con- 
stitute the solution to the basic equations of § [2] By looking 
at the solution, we can identify cases where a cavity is kept 
empty in steady-state or when the disk overflows. 

The cavity refills if the inner edge of the disk outside the 
secondary, n, falls within the Hill radius of the secondary. 
The viscous torque in the near zone, either T° es or T" eu , is a 
monotonically decreasing function of n. Thus, if the disk is 
truncated, then the viscous torque at any radius in the near 
zone is decreased relative to its value for an overflowing disk, 
n = r s + 7"H (Eq. l82|) . This shows that the state of the disk is 
uniquely determined by the smallest torque barrieiF"!. Tb c , or 
equivalently, the smallest dimensionless angular momentum 
flux: 

fc s = min{fc s m ° u , fc s msu , fcr gs } ■ (134) 

given by Eqs. (|116H117fl and (|130H131[) . We therefore must 
distinguish four different possible cases of migration and disk 
behavior. First, fc s = fc™ gu corresponds to the standard case 
with a wide gap, with the secondary exhibiting Type-II mi- 
gration. If fc s = k™ ss , then the gap edge is located within a 
scale-height in the near zone so that the torque cutoff limits 
the tidal torques, but they can nevertheless support a gap 
against viscosity as the secondary migrates inward. How- 
ever, if fc s = k™° B , then the saturated tidal torque becomes 
smaller than the viscous torque all the way to the Hill ra- 
dius, and the cavity refills. Finally, if fc s = fc™ ou , then the 
disk reaches the Hill radius and overflows already while the 
torques in the near zone are still unsaturated. 

The migration rate of the secondary is proportional to 
fc s , (Eq.HHJ, i.e. 

|w Br | = 2k s Mr s /m s = min{«° u , C, «fr} ■ (135) 

given by Eqs. (|120H12ip and (|132H133|I . Note that wf r u and 
wfr are given by practically the same formula, up to an 
order-of-unity factor of A. This correspon ds to the case 
of se condary-dominated Type-II migration (|Sver fc Clarke! 
Il995l ). The migration rate in the overflowing case is «™ or 
for unsaturated or saturated torques. Here the disk is 
still strongly perturbed, but gas inflow across the orbit lim- 
its the efficiency of migration. The disk structure in this new 
regime is intermediate between a disk with an empty gap 
(normally associated with Type-II migration) and a weakly 
perturbed disk (Type-I migration). Although the migration 
speed in this regime is slower than either in standard Type-II 
or Type-I migration, we refer to this regime as "Type-1.5". 

The gap can also close due to three dimensional overflow 
for large secondary masses if the tidal heating is substantial 
to make the disk puff up. Eq. (|56p shows that this happens 
(H > r) in a radiation pressure dominated disk if 

S thick r^j-r 4/15 . 1/3 „ r-2/b , -7/6 1/3 /,,m 

r s $ T s = 57A/. a_ 1 mj_ x M 7 ' A ' qj 3 (136) 

or equivalently if 

q > qthick = 5.4 x 1CT 3 aZ^'mZlM^ 5 A 7/2 r s 3 2 . (137) 

In this case, we do not derive the geometrically thick over- 
flowing disk or the migration rate. 

20 recall that T^ c = T v at the interface between the near and 
middle zone 



We discuss the gap opening and closing conditions in 
more detail in Paper II, where we contrast them explic- 
itly with the standard expressions widely used in the litera- 
ture, and also compare Type-1.5 migration to the previously 
known Type-I and II cases. 



5 RESULTS AND DISCUSSION 

We have derived analytical solutions to the disk model in 
different radial regions, where either the tidal, the viscous 
torques, or the angular momentum flux is negligible relative 
to the other two terms in Eq. In particular, we have 

identified the far zones, either well inside or outside the sec- 
ondary's orbit, where the effects of the secondary are negli- 
gible, the exterior middle zone, where the disk structure is 
greatly modified but where the tidal torque and heating are 
locally negligible, and the near zones just inside or outside 
the secondary, where the tidal effects dominatd 2 *!. We distin- 
guish two cases in the middle zone, depending on whether 
the disk has a gap (i.e. the disk is truncated well outside 
the Hill radius) or if the disk is overflowing across the sec- 
ondary's orbit. We also distinguish two cases in the near 
zone outside the secondary's orbit, depending on whether 
the tidal torque from the binary is saturated (i.e. whether 
the location of the cavity edge falls within a scale height 
H from the secondary). Furthermore, we have investigated 
asymptotic results for gas and radiation pressure dominated 
cases. 

Distinguishing all of the above cases allowed us to adopt 
separate approximations, each valid in the corresponding 
regime, and to derive analytical results to the perturbed 
accretion disk interacting with the secondary. We further 
used this to estimate the migration speed for the secondary. 

In this section, we collect all of the resulting analytical 
solutions for the most important disk parameters in the var- 
ious zones, and present them in a form suitable for easy use. 
We refer the reader to Paper II for physical interpretations, 
and discussions on possible implications of our results for 
real binary systems. 

5.1 Disk model 

Our results in this paper apply to geometrically thin, opti- 
cally thick accretion disks, and describe vertically and az- 
imuthally averaged properties. All physical parameters can 
be written as 

X(r, r „ p) = C a c _\ rh c _\ r? f% q% r% *(r, r„ p) 

(138) 

where X denotes any of {E, T c , H, v r , F, T„}; ri and r S 2 
denote the radius from the primary and the orbital ra- 
dius of the secondary in units 100 M., respectively; and 
$(r, r s ,p) denotes an extra function of the parameters p = 
(a,rh, f, M,q). Note that (a,m,M) are the usual parame- 
ters of a standard solitary accretion disk; q and / represent 
the mass ratio and the normalization of the azimuthally 
averaged tidal torque (see Eq. [5|. The _jv index denotes 



21 We have further restricted the near zones to lie outside the 
Hill radius of the secondary. 
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Table 2. Pre-factors and exponents in the analytical disk model in different zones, C and Ci in Eq. I|138|l . The third column is in cgs 
units except where marked by * (where it is in units of GM, /c 2 ). Columns 4-10 are exponents, the last column is the extra multiplicative 
function (see text). The last two block of parameters show the migration rate of the secondary and other useful parameters. 
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C [a-i] [m_i] [My] [/_ 2 ] [ g _ 3 ] [r s2 ] g 

Middle/Far w. gap sat. r™ gs /r s 540 1 -3/4 -3/2 574 -7/4 (1 + <5£ 1 ° g s )~ 11/8 

Middle/Far w. gap uns. rf eu /r B 540 1 -3/4 -3/2 5/4 -7/4 (1 + <S£i C g u ) _11 / s 

Near ext./middle w. gap sat. 5£f| 0.130 -4/37 4/37 -3/37 9/37 22/37 -33/74 

Near ext./middle w. gap uns. Sgg 0.55 -4/37 -6/37 -3/37 9/37 22/37 -13/74 (1 + i^")- 35 / 148 

Near ext./middle w/o gap sat. 6™ 0.083 -4/17 1/17 5/17 10/17 -5/34 |W| 4/17 

Near ext./middle w/o gap uns. 6™£ 0.35 -4/7 -2/7 1/7 5/14 20/21 -1/14 (r ; /r s )- 115 / 126 

Far/Near int. S { ni -0.1 -4/17 -2/17 1/17 5/17 10/17 -1/34 (1 + <5f ,)-84l/7l4 

Table 3. Transition radii between different zones relative to the secondary orbital radius, r B . Here is the radius at the interface 
between zone a and b, <5rjj = r Jj — r B , and S b = 5r£ /r B . Different columns show the constant prcfactor and exponents in Eq. d 1 38 D as in 
Tabled 



normalizing with 10 _JV , e.g. g_3 = g/10~' ! . The C prefac- 
tor, c, exponents, and the $ function in Eq. ()138p are given 
in Table [2] in the different zones and cases for /3-disks (i.e. 
v oc Pgas). For /3-disks, many of the physical parameters, in- 
cluding the surface density, E, and the central temperature, 
T c , are independent of whether gas or radiation pressure 
dominates. This, however, is not true for the scale-height, 
where we quote results in both regimes, labeled with a "gas" 
or "rad" subscript. 

The last column of Table [2] shows $(r, r s ,p) in 
Eq. (I138[) . Here we introduced the following notation: 



8r 
A, 



r — r s , Sr\ 



Sfl = (fi 8 - Q)/fi s = 1 - (r/r s y 3/2 , 
Silo = (O s - fio)/U» = 1 - (n/r s )" 3/2 . 



1/2 



(139) 
(140) 
(141) 
(142) 

(143) 



/2 fr/r s 



5/2 , 



-7/6 



(l-^ 2 )" 1/6 (a; -l) 



2 1 / 6 2 (n/r,)- 115 / 72 
3 1 / 6 5 (o>i/r H ) 5 / 2 



r .s 115/72 /A; 



(?) 



-10/3 



5/2' 



(144) 



\ 

J 5SI 



-21/26 



(1-X) 



-6/13 



d.r 



= |s(d~ft; 5/26, 7/13) - ~-B(d"fi ; 5/26, 7/13) 

«0.41nf", 

on 

W = W-i(-a) ~ ln(a) - ln(- ln(a)) , 
-us/ 72 / Sr . \ -a/ 2 



(145) 
(146) 

(147) 



Here B(a;; a, b) is the incomplete Beta function, and W is 
the Lambert-W function, the branch defined on the nega- 
tive real axis, evaluated at —a given in the last row of Ta- 
ble H Typically, 1 < |W(-a)| < 10. The approximations 
shown for ip and W are better than 20%. Here, nsco is the 
innermost stable circular orbit near the SMBH (i.e. 6M. 
(1M.) for a non-spinning (maximally spinning) SMBH)), 
th = (g/3) 1//3 r s is the Hill radius, n marks the radius at 
which the viscous torque becomes very small outside the 
secondary's orbit for which we use 



and A is the dimensionless gap or truncation radius scale 
\ As if k B = kT SB , ( ' 

wherH 



k s — {k s 



(150) 



Sri — th if k s = fc™ " or k s = fc s 



(148) 



is the dimensionless angular momentum flux or brightening 
factor. 

We collect the formulae for the transition radii separat- 
ing different radial zones and physical regimes in Table 
We label = (r„ — r s )/r s , where r b a marks the radius of 
the interface between zone a and b. For example, S^g is the 
transition between the middle zone and the torque-saturated 
exterior near zone if there is a wide gap, which also sets the 
truncation radius scale A s in Eq. (|140[) . 

The state of the disk and the migrate speed of the binary 
is directly set by the dimensionless angular momentum flux 
k B . If k B =[k™ ou ]kT OB , then the disk is in the [un] saturated 
overflowing state, whereas if k B =[fc™ su ]fc" lss then it is in the 
[unsaturated state with a wide gap. These four possibilities 
are indicated in Table [2] The appropriate choice of k B also 
determines which case in Tableware to be used in the middle 
zone for the other parameters ('mou', 'mos', or 'mg'), and 
in the near exterior zone ('neu' or 'nes'). 

Thus, the "phase space" of solutions consists of four re- 
gional- The transition between two different solutions, with 
or without wide gaps, corresponds to the parameters for 
which k™ os = ££ lgs or k™ ou = fc s ngu - Our hypothesis is that 
this must represent a physical transition in the disk+binary 
system, as the secondary migrates inward from large radius. 
Initially, a central cavity is created, and the outer edge of the 
cavity lies far from the secondary's orbit. However, as the 
secondary migrates inward, the distance between the cavity 
edge and the secondary shrinks (at least when measured in 
units of the Hill radius of the secondary). This may happen 
both because the viscosity increases as the pressure grows 
during pile-up, and also because the tidal torque decreases 
with increasing scaleheight due to the torque cutoff. The 
cavity finally closes once the cavity wall nudges inside the 
Hill radius. The dimensionless angular momentum flux or 
brightening factor, k B , is largest when this transition occurs. 



22 If either kT SU , k™ ss , fcf ou , or fcf os is less than one, then the 
disk does not have a middle zone by definition, and the disk is 
not strongly perturbed. 

23 Here we assume that GW emission is negligible, and the disk 
drives the binary. See Paper II for further discussion. 
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Figure 2. Comparison of the numerical solution (thick yellow solid lines) with the asymptotic analytical approximations in the various 
zones (dotted, dashed, and dash-dotted lines, labeled as in Table [2)| for the surface density of the disk around a 10 5 Mq primary. The 
mass ratio and binary separation are (q,r B ) = (0.01, 500 M) and (0.1, 1000 M) on the left and right panels, respectively. In both cases 
the disk is in the overflowing steady-state. The tidal torque is unsaturated in the near zone outside the Hill radius on the left panel, but 
it is saturated on the right panel. The disk structure is significantly modified in both cases in an extended region around the secondary. 



These gap opening/closing co nditions are quite differ- 
ent from those in the literature (e.g. I Ward! 19861 ) . which state 
that the gap closes if the disk extends into the region closer 
than either the local scaleheight or the Hill radius from the 
secondary. Note that our gap closing condition combines 
statements on the scaleheight and the Hill radius, but since 
the scaleheight and viscosity vary significantly near the sec- 
ondary in the strongly perturbed case with a large pile-up, 
they depend on the actual perturbed profiles rather than 
the averaged quantities describing accretion disks around 
a solitary object. We discuss gap opening/closing, and its 
physical implications, in more detail in Paper II. 

We have verified that the analytical approximate solu- 
tions to the disk model match the numerical solutions typi- 
cally to within tens of percent for a wide range of disk and 
binary parameters when the disk is strongly perturbed. Fig- 
ure [5] shows two examples when the disk is in the overflow- 
ing state with unsaturated (left panel) and saturated (right 
panel) tidal torques near the secondary. Different regions are 
indicated with the abbreviations used in Tabled For further 
examples and other physical quantities, see Paper II. 

5.2 Brightening factor 

As stated in the previous sections, the k s parameter sets 
the angular momentum flux and the brightness of the disk 
in the middle zone relative to the unperturbed value, and 
determines the state of the disk. 

Remarkably, the disk parameters can differ dramati- 
cally from the unperturbed values not only in the near zone, 
where the tidal effects dominate, but also in the middle zone, 
where tidal effects are already negligible. The tidal effects 
of the secondary are short-range, but the corresponding ef- 
fect is communicated to distant regions by setting an ef- 
fective boundary condition. The size of this region can be 
r f mg /r s ~ 540a-irhZi^ 4 M^ 3 ^ 2 r~2^ 4 qt_{j times larger than 
the secondary orbital radius (see Table |3J. This long-range 
behavior is co nfirmed in our full n umer ical solutions and is 
also present in lLodato et all |2009h and lChang et al.1 (|201Ch 



for a circumbinary disk with a cavity. The same k s param- 
eter also sets the increase in the scale-height, as well as the 
migration speed of the secondary. We have derived the ana- 
lytical formulae for k s in the four relevant regimes summa- 
rized in the previous section. Due to its extended radial size, 
the integrated luminosity of the middle zone is often larger 
than that of the near zone in the torque-unsaturated case, 
overflowing state. However, the near zone may be brighter 
than the middle zone, when the system is in the torque- 
saturated, overflowing state (see Paper II for further discus- 
sions) . 

5.3 Migration rate 

Finally, we have derived the migration speed of the sec- 
ondary; approximate analytical formulae for our results are 
listed in Table [2] We find that in the regimes when a cavity 
forms, the migration in our solution is slower than the stan- 
dard steady-state second ary-dominated Type-II migration 
rate (|Sver fc Clarkelll995l ) by a factor A 11/16 ~ 2. The Type- 
II rate may be further reduced in non-steady state models 
if the accretion rate or th e gas mass is limited (|lvanov et al.l 
1 19991 ; lLodato et al.1 [20091 ). However, in the new "Type- 1.5" 
regime we identify, with a partial pile-up and overflow, the 
migration is even slower. 

More generally, the migration speed for a strongly per- 
turbed disk follows the minimum of the three possible solu- 
tions: 

|« sr | = mm(|u sr ,1.5s|, |Usr,1.5u|, |«sr,Il|) ■ (151) 

In general, Type-1.5 migration is more rapid at larger 
r s , in contrast with the Type-II rate, which is nearly con- 
stant, and the GW inspiral rate, which strongly decreases 
with r s . Therefore, as a real system evolves, it will first 
transition from an initial Type-II migration at large radii 
to Type-1.5 migration at smaller radii, before finally being 
driven by GWs at still smaller separations. 

Regarding the dependence on the primary mass, the 
Type-II is nearly constant, while the Type-1.5 speed in- 



© 0000 RAS, MNRAS 000, 000-000 



Gas pile-up, gap overflow, and Type 1.5 migration: general theory 19 



1000 



o 



T3 



10 



0.1 



Type II 

numerical solution 



analytical prediction 
__, i . . . . i i ■_ 



1000M. 

© 



M. = 10 5 M G 




0.001 



0.01 



0.1 



Figure 3. The migration speed of the secondary for different 
mass ratios for M. = 10 5 Mq at r B = 1000M. (or t orb = lday). 
The steady-state numerical solutions (blue) are well represented 
by our analytical formula for Type- 1.5 migration (magenta), but 
not by those for Type-I or II (black dotted). The dashed and 
solid magenta lines show the unsaturated and saturated Type-1.5 
cases, respectively. The discrepancy between the numerical and 
analytical solutions become more significant at q < 0.002 due to 
the fact that the disk is not strongly perturbed there. 



creases with M # . This implies that at any fixed orbital 
separation, Type-1.5 migration is relevant for lower masses 
(roughly those in the range expected to be detectable by a 
space-based gravitational wave mission such as LISA, and 
Type-II is relevant for higher-mass binaries (in the sensitiv- 
ity range of Pulsar Timing Arrays; see Paper II). Type-1.5 
migration may also be important for stellar mass binaries in 
proto-stellar disks or Jupiter-mass planets around M-t ype 
dwarf stars of mass (0.1-1) Mq (see lJohnson et aT1l2012l . for 
a recent discovery of such a system). Finally, we note that 
Type-1.5 (Type-II) migration operates for larger (smaller) 
accretion rates, if fixing all other parameters. 

All of the above conclusions are based on the analyti- 
cal solutions we obtained; however, we have verified, by nu- 
merically solving the equations presented in § [5] that our 
solutions are accurate to within tens of percent for a broad 
range of parameters. In particular, in Figure [3] we show the 
the analytical approximation of the migration rates for the 
case of r s = 1000 M. for M. = 10 5 M Q , together with the 
rates obtained from a numerical solution. In this disk 
is in the overflowing state for all values of q > 10~ 3 shown, 
and the migration rate is significantly different from both 
Type-I and II. As the figure shows, the analytical Type-1.5 
formulas give a good approximation over a wide range of q 
for the migration speed. 



5.4 Caveats 

Our findings are subject to many possible caveats. 

• We assumed a radiatively efficient disk model in which 
the effective viscosity is proportional to the gas pressure in 
the disk with a constant a coefficient even in the radiation 
pressure dominated regime. Future studies should investi- 
gate alternative models in which the vis cosity is proportional 
to the total gas+radiation pressure jShakura fc Sunvaevl 



Il973l ). or where the viscosit y is generated by magneto- 
rotat i onal instability (MRI, seetTurner e t al. 2003; Sh i et all 
|2012| ; iGiacomazzo et al.ll2012l ; iNoble et al.ll2012l for simula- 
tions of circumbinary disks leading to an "antigap"). 

• We assumed steady-state models where the accretion 
rate is constant over radius. This is expected to be valid as 
long as the gas inflow is much faster than the migration rate 
of the secondary over a large range of radii, if the total gas 
supply is not limited, and if the accretion rate is set at the 
inner or outer boundary (i.e. for the new Type-1.5 migration 
regime we focus on here). However, this assumption may be 
violated fo r tidally truncate d circumbinary disks (Type-II 
migration. Ilvanov et al.lll999l ) or for models where the gas 
supp l y rate is limit ed at the outer boundary |Lodato et al.l 
l2009l ; lRafikovll2012l ). 

• We assumed unequal-mass binaries, averaged over the 
azimuthal angle, assumed that the density waves gener- 
ated by secondary are dissipated locally, and that the 
radial tidal torque pr o file f ollows the formula given by 
lArmitage fc Natarai an (2002). This assumes that the tidal 
torque saturates near the secondary at a radial distance 
closer than the scaleheight {d r T^ oc H~ 4 ). We also assumed 
that the gas entering a distance comparable to the Hill ra- 
dius can flow freely across the secondary's orbit. However, 
accretion onto the secondary may affect the Type-1.5 mi- 
gration rate. Farther away from the secondary, the assumed 
torque density has a steep cutoff (<9 r Td oc \r — r s ]~ 4 ); extrap- 
olating beyond r > 2r s might be inaccurate. The |r — r s |~ 4 
scaling may also be inaccurate in the local nonlinearly per- 
turb ed regime especially for compara b le mass-ratio bina- 
ries dMacFadven fc Milosavlievicl 120081 ; iRoedig et all |2012| ; 
IPetrovich fc Rafikovl |2012| ). These issues should be inves- 
tigated using simulations, which could also address com- 
parable mass binaries w here the disk may be significantly 
non-axisymmetric (e.g. iMacFadven fc Milosavlievicl 120081 ; 
ICuadra et ail 120091 ) and where the accretion of the sec- 
ondary is non-neglig ible (|Lubow et al.lll999l ). 

• We neglected non-axisymmetric inflow into the sec- 
ondary's orbit or onto the secondary if a cavity is formed. 
Inflow across the gap or accretion of the secondary reduces 
the amount of pile up outside the secondary, reduces the 
Type-II migration rate, and could affect the gap closing tran- 
sition between the continuously overflowing solutions and 
the cases with a gap. We have also neglected the corotation 
torques in the overflowing case. 

• We found that the enhanced pressure dominates over 
the increase in surface density outside the secondary's or- 
bit, which makes the overflowing disk stable against gravita- 
tional fragmentation (see Paper II). However, the steep pres- 
sure gradient in the near zone around a massive secondary 
may lead to global non-axisymmetric dynamical ins tabili- 
ties (|Papaloizou fc Pringlelll985l ; lGoldreich et al.lll986t ). The 
corresponding enhancement of the effective viscosity and an- 
gular momentum transport in the disk might reduce the 
pile-up outside the secondary's orbit and further reduce the 
Type-1.5 migration rate. A detailed stability analysis and 
an investigation of its implications goes beyond the scope of 
this paper. 

• We restricted our attention to circular binaries. How- 
ever, binary eccentricity may be excited in cases with a cav- 
ity, when the masses of the two compact objects are com- 
parable and the gap edge itself becomes significantly non- 
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axisymmetric (Artymowicz _ 19921 ; lArmitage fc Nataraianl 
120051 ; iRoedig et al.ll201ll , |2012h . We note that such non- 
axisymmetries excited in disks with a gap diminish once the 
mass ratio is q < 0.1 (D'Orazio et al. 2012, in preparation). 
Nevertheless, it remains to be seen if the binary develops sig- 
nificant eccentricities in the unequal-mass, overflowing state, 
with a significant pile-up. 

• We assumed that the binary is sufficiently widely sep- 
arated that gravitational wave emission is negligible. For 
a complete picture, future studies should investigate the 



gravitational wave driven regime dTanaka & Menoul 


201C ; 


Chang et all |2010|; lYunes et al.1 l201ll; iKocsis et alJ 


2011; 


Farris et al. 201 ll; Tanaka et al. 2012: 


Bode et alJ 


2012; 


Baruteau et al.l 20121; Giacomazzo et al.l |2012|; Noble et al. 


2012D. 



6 CONCLUSIONS 

In summary, we have presented new analytical solutions 
to disk properties and migration rates, obtained from self- 
consistent solutions of a coupled binary-disk system. The 
evolution equations are solved analytically in the strongly 
perturbed limit, including the angular momentum exchange 
between the disk and the binary and the modifications to 
the density, scale- height, and viscosity self-consistently, in- 
cluding viscous and tidal heating, diffusion limited cooling, 
radiation pressure, and the orbital decay of the binary. 

In addition to recovering solutions with a central cavity, 
similar to previous "Type-II migration" scenarios, we have 
identified a distinct new regime, applicable at smaller sepa- 
rations and masses, larger accretion rates, and mass ratios 
in the range 10 -3 < q < 0.1. For these systems, gas piles up 
outside the binary's orbit, but rather than creating a cav- 
ity, it continuously overflows as in a porous dam. The disk 
properties are intermediate between those in an unperturbed 
disk and a disk with a wide gap. The migration rate of the 
secondary in this "Type 1.5" regime is typically slower than 
both Type-I and Type-II rates. 

In this paper, we have presented simple analytical for- 
mulae that comprehensively describe binary systems with 
different parameters, in various stages of evolution. The an- 
alytical results provide simple scaling relations, which may 
be useful to scale and interpret the results of numerical sim- 
ulations to different disk or binary parameters. It allows us 
to map out the effects of varying a, m, or the binary pa- 
rameters, over a wide range from planetary disks to active 
galactic nuclei around SMBH binaries. 

We discuss the applicability of the new Type- 1.5 regime 
and its physical implications for specific systems, as well 
as possible observ able signatures in a companion paper 
(|Kocsis et al.ll2012h . 
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APPENDIX A: VISCOUSLY AND TIDALLY 
HEATED DISKS 

Here we provide the details of the algebraic manipulations 
that give the disk model for fixed F and D v for either a or 
/3-disks (i.e. b = or 1). 

For fixed /3, we can reduce the problem to 3 equations 
and 3 unknowns E, T, v. From Eqs. |15J), and (fl8|) . 



8 k 2 p b F 2 , 

= 9Q 2 = a Q^Jl~PY' 

3 K E 



Solve this for E and T, 



E = 



8c J 



■Q 



9ctK 2 

t: = -— n 



Rb pi 



Sana" /3 b F 

Finally f3 is given by 

P _ Pgas _ 3pkT c 3k E 

1 — ft prad aiirripTc afim p 2HTj? 

where E, T c , and H are to be substituted from Eqs. (I A3I) . 
(A4|l . and (16]). From this 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 



P 



(l/2)+(6-l)/10 



1-P 



c[fc/( M m p )] 2 / 5 ^ /w Dl /W 



Finally Eqs. (|A3[) and (|A4[I can be simplified by eliminating 
1-/3 using Eq. CM}, which leads to Eq. (|21H22I) . 
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